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Abstract  / 

Progress  is  reported  for  the  contract  period  covering  the  past  18  months. 
Substantial  advances  in  the  analysis  of  the  large-amplitude  unsteady  aerodynamics 
of  wings  have  been  achieved  and  are  summarized  in  the  attached  papers.  'Special 
focus  has  been  placed  on  developing  techniques  for  describing  the  exact  non-linear 
convection,  deformation,  and  roll-up  of  the  vorticity  wakes  above  and  behind 
active  wing  surfaces.  Using  a  combination  of  analytic  and  computer-interactive 
methods,  new  insights  and  more  exact  aerodynamic  , performance  results  have  been 
generated.  The  effects  of  dissipative  internal  "cores"  in  the  overall  3D  wake 
structures  above  deltas  have  been  analyzed  in  the  slenderj-wing  limit.  New  paths 
for  further  computer-assisted  advances,  including  improved  understanding  of 
wing-wake  interaction  in  severe  maneuver,  vortex  core  formation  and  break-up, 
and  aerodynamic  history  effects  have  been  charted.  '>■  e  v 


Statement  of  Progress 

The  papers  attached  summarize  as  carefully  as  we  are  able  the  sequential 
aspects  of  the  progress  made  during  the  above  contract  period.  In  addition  to 
these  three  papers,  two  key  Master  of  Science  thesis  reports  (see  attached  list) 
have  proved  to  be  major  contributors  to  our  understanding  of  wing  aerodynamics 
involving  both  large-scale  unsteady  motion  and  viscous-related  phenomena  such 
as  leading-edge  separation,  bubble  formation  and  reattachment  and  vortex  core 
evolution. 

In  paper  1  "Nonlinear  Aerodynamics  of  Two-Dimensional  Airfoils  in  Severe 
Maneuver"  (Scott /McCune) ,  advanced  computer-interactive  techniques  for 
studying  large-amplitude  unsteady  wake  motion  behind  host  airfoils  is 
described.  The  nonlinear  generalization  of  the  classical  Wagner  equation 
is  included. 
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In  paper  2  "Interactive  Aerodynamics  of  Wings  in  Severe  Maneuver?'  (McCune, 
Colorado  Springs),  analysis  of  exact  nonlinear  unsteady  wake  convection 
over  slender  delta  wings  in  severe  maneuver  is  described.  Application  of 
paper  1  to  that  problem  is  detailed. 

In  paper  3  "Slender  Wing  .Theory  Including  Regions  of  Embedded  Total  Pressure 
Loss" _ (McCune ,  Tavares,  Lee,  Weissbein) ,  self-consistent  wake  development  over 
slender  wings  including  the  development  of  dissipative,  low-total-pressure 
core  regions  is  described.  Viscous  related  criteria  for  the  sizing  and 
evolution  of  core  size  shape  and  intensity  are  discussed.  The  importance 
of  the  "history"  effect  of  wake  development  on  wing  performance  is  pointed 
out . 
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Abstract 


Yc  (x.t) 


This  paper  presents  a  nonlinear  theory  of  forces  and 
moment  acting  on  a  two-dimensional  airfoil  in  unsteady  uf  8) 
potential  flow.  Results  are  obtained  for  cases  of  both 
large  and  small  amplitude  motion.  The  analysis,  which  uw  ( 81 
is  based  on  an  extension  of  Wagners  integral  equation 
to  the  nonlinear  regime,  takes  full  advantage  of  the 
trailing  wake’s  tendency  to  deform  under  local  veloci-  uw(6i 
ties.  Interactive  computational  results  are  presented 
that  show  examples  of  wake-induced  lift  and  moment 
augmentation  on  the  order  of  20  percent  of  quasi-static  s 
values.  The  expandability  and  flexibility  of  the  present 
computational  method  are  noted,  as  well  as  the  relative 
speed  with  which  solutions  are  obtained.  smaz  ( t) 
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hit) 


Nomenclature 
Angle  of  attack 

Height  of  airfoil  midchord  above 
reference  axis 


Yu/(».  t) 
UdB» 


Position  of  the  airfoil  or  camber  line  in 
physical  (airfoil)  plane 

Total  tangential  velocity  on  the  airfoil 

Wake-induced  tangential  velocity  on 
the  airfoil 

Wake-induced  velocity  normal  to  the 
airfoil 

Distance  along  the  wake  from  the 
trailing  edge 

The  location  of  the  end  of  the  wake 
farthest  from  the  trailing  edge 

Strength  of  vorticity  in  the  wake 

Average  tangential  (parallel)  velocity 
component  on  the  airfoil: 

{u!  Bi)  =  +  ufB  +  it))/2 
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8  ,  T 
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Radial  and  tangential  components  of  Hu(B) 
velocity  in  the  transformed  (circle) 
plane 


Radius  of  airfoil  in  transformed  (circle) 
plane 

Angular  running  variables  in  trans¬ 
formed  (circle)  plane 

Circulation  strength 

Cartesian  coordinates  fixed  on  global 
reference  axes  in  the  physical  (airfoil) 

plane 

Cartesian  coordinates  parallel  and 
normal  to  the  airfoil 
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r0  (t> 

Cc  t  > 


The  difference  in  velocity  components 
on  the  upper  and  lower  surfaces  of  the 
airfoil;  equal  to  the  strength  of  the 
vorticity  on  the  airfoil: 

Sui B)  =  id B)  -  id B  +  n)  -  y(B) 

Free  stream  velocity  and  density, 
respectively 

Potential  function 

Lift  and  moment  on  the  airfoil 

Chord  length 

Quasi-steady  bound  circulation 
Discrete  vortex  location  in  the  wake 
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x  i  P  Location  of  the  trailing  edge 

Introduction 

Early  studies  of  the  linearized  response  of  the  forces 
and  moment  acting  on  an  airfoil  in  unsteady  flow, 
including  the  effects  of  its  own  unsteady  wake,  have 
been  described  by  such  authors  as  Theodorsen.'  von 
Karman.2  Kussner.J  Sears. <  and  Wagner.5  Limiting 


assumptions  in  the  classical  two-dimensional  analysis 
are  that  the  airflow  be  idealized  as  a  potential  flow,  that 
the  wake  be  approximated  as  trailed  out  flat  behind  the 
airfoil,  and  that  the  amplitude  of  airfoil  motion  be 
correspondingly  small.  However,  modern  requirements 
for  aircraft  maneuverability  now  dictate  that  problems 
associated  with  fully  nonlinear,  large-amplitude, 
unsteady  flows  be  addressed.  Thus,  the  theoretical  basis 
for  the  linearized  response  must  be  reformulated  and 
extended  to  allow  the  inclusion  of  nonlinear  effects, 
including  wake  deformation.  The  wake  mus£  be  allowed 
to  convect  freely,  so  that  an  effective  "history"  of  the 
airfoil  motion  is  present  in  the  wake  elements.  An 
extension  of  Wagner's  integral  equation. 5  is  necessary 
to  relate,  in  a  time  accurate  sense,  the  strengths  and 
positions  of  each  of  these  wake  elements  to  the 
circulation  on  the  airfoil. 

With  the  interactive  capabilities  of  a  VAX  11-750 


In  equation  (4),  Yc  ix.t)  denotes  the  position  of  the  plate 
or,  more  generally,  the  camber  line  of  a  given  airfoil. 
Given  this  formulation,  we  note  that  the  substantial 
derivative  D^/Dt  must  vanish  at  y=0. 10  Therefore,  if 
the  equation  for  the  position  of  the  airfoil  is  defined,  the 
velocity  on  the  airfoil  at  a  given  time  t  and  position  x 
may  be  computed.  The  present  formulation  of  Y^x.t) 
involves  the  uncoupled  use  of  two  parameters:  Aft),  the 
height  of  the  midchord  above  the  reference  axis,  and 
<xtti.  the  angle  of  attack  relative  to  the  free  stream.  An 
input  or  "forcing"  function  is  defined  by  setting  the 
variation  of  these  two  parameters  with  time.  The  input 
function  then  yields  surface  velocity  components 
through  the  use  of  the  substantial  derivative.^ 

Conjugate  function  theory  requires  that  the  veloci¬ 
ties  impinging  upon  the  airfoil  due  to  the  presence  of 
wake  vorticity  be  "removed”  from  the  problem  before 
proceeding.  This  process  is  illustrated  in  Figures  1 


computer  and  solving  for  only  one  unknown  per  time 
step,  it  is  possible  to  convect  the  wake  very  quickly,  and 
this  enables  the  user  of  the  program  to  see  the  wake 
unfold  and  flow  downstream  as  the  program  runs. 
Output  of  force  and  moment  histories  is  generated  on¬ 
line.  enabling  a  quick  test-and-evaluation  turnaround 
time  of  only  a  few  minutes  per  case.  The  following 
sections  outline  the  derivation  and  implementation  of 
this  interactive  nonlinear  solution  and  present  results 
for  typical  input  histories. 


through  3. 


Derivation 

Laplace's  equation  in  the  vicinity  of  an  oscillating 
airfoil  in  a  two-dimensional  flowfield  can  be  solved 
explicitly  by  means  of  the  classical  Joukowski 
transformation  of  the  region  to  complex  space  .1'  By 
using  the  theory  of  conjugate  functions,  and  following 
the  explicit  derivation  of  reference  8,  we  can  describe 
the  radial  and  tangential  components  of  velocity  that 
exist  on  the  airfoil  in  the  transformed  plane  as  follows: 


Figure  1.  Apparent  wake-induced  velocity  on  an  airfoil. 


Figure  2.  Velocity  at  airfoil  surface  due  to  additional 
bound  vorticity. 


v*'a-ll‘s-kh  tvaf'ro'(Lr)d“2-^  ,2) 

To  generate  the  appropriate  circulation  on  the  airfoil 
and  in  the  wake,  the  Kutta  value  of  circulation  is 
applied,  and  (2)  can  be  rewritten: 


I  X  sin  3  +  sin  i  Is-  <<.,  ( 

u,  !i  •  -  -  —  If  v  (a.  1 1  - - -  -  -  (I  i 

2n  ]<i  ’  cos  3  -  cos  c  ,in  i 


Boundary  conditions  at  the  surface  of  the  airfoil  are 
prescribed  in  the  manner  of  reference  9  by  the  use  of  an 
implicit  variable  formulation.  For  purposes  of 
simplicity,  the  airfoil  in  the  present  study  has  been 
assumed  to  be  a  fiat  plate  whose  instantaneous  position 
relative  to  a  global  reference  system  can  be  described  by 
a  function  q  as  follows: 


1  x.  v.  t  -  v  -  Y  x,  t  - 


(4) 


x 


Figure  3.  N'et  flow  with  normal  velocity  cancelled  at 
airfoil  surface. 

Figure  1  shows  the  apparent  wake-induced  velocity 
impinging  on  the  plate.  In  order  to  counter  the  normal 
velocity  a  "fictional"  problem  must  be  created  (Figure  2) 
in  which  additional  bound  vorticity  is  produced  on  the 
surface  of  the  airfoil.  When  these  two  solutions  are 
added  (Figure  3),  the  no-flow-through  boundary  con 
dition  of  the  net  flow  on  the  plate  is  achieved.  Conjugate 


function  theory  must  create  the  "fictional”  flow  by 
measuring  the  wake-induced  upwash  at  each  station 
and  then  inducing  the  bound  vorticity  on  the  plate  to 
counter  the  upwash.  Thus  we  may  write  that  the  total 
velocity  component  parallel  to  the  plate,  u( &),  is 
composed  of  the  addition  of  the  wake-induced 
component,  uJB),  and  the  conjugate  function  denoted 
by  equation  (3). 
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uj  B)  may  be  derived  by  the  use  of  the  Biot-Savart 
law  in  complex  form: 


Calculation  of  loading,  lift,  and  moment  proceeds 
along  classical  lines  by  means  of  the  unsteady  Bernoulli 
equation: 


-*p=p. 
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(9) 


Reduction  of  the  resulting  expressions  requires  a 
judicious  use  of  the  Glauert  Integral  and  ultimately 
leads  to  the  following  final  equations:* 
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(6) 


Since  the  circulation  on  the  airfoil  can  be  described 
as  the  difference  between  the  velocities  uf8)  and 
u/B  +  n),  the  loading,  lift,  and  moment  can  be  computed 
if  the  velocity  u(  8)  can  be  found  for  values  of  8  between  0 
and  2ir.  One  very  important  outcome  of  this  analysis  is 
that,  for  the  general  nonlinear  case,  the  following 
expression  is  always  true: 


(  u  ip> )  =  U m  cosa  (  n  *■  u^tp) 


(7) 
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where  (u(BJ>  is  taken  to  be  the  average  of  the  velocities 
ti(  8)  and  u  ( 8  +  v).  This  result  proves  that  the  only  in¬ 
fluencing  factors  on  the  average  platewise  velocity  are 
the  free  stream,  the  angle  of  attack,  and  the  location  of 
the  wake's  vortex  elements.  The  linearized  theory  is 
recovered  if  the  wake  is  trailed  out  fiat  behind  the  airfoil 
and  the  angle  of  attack  is  small.  Equation  (7)  then 
reduces  to  the  familiar  expression  iu(B))  =  U.  as  re¬ 
quired.2 


The  expression  for  the  vorticity  on  the  plate  at  a 
position  8  is  more  complicated,  but  still  conveys 
meaning: 


-  y  *0)  -  u  (p)  —  u  ifi  +  n)  =  2U *sin  a  W )  —  2(J acosa  u 
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where  uu,U)  is  measured  at  a  station  £  on  the  airfoil  and 
is  derived  by  means  of  equation  (6).  We  notice  that  the 
first  two  terms  on  the  right-hand  side  of  (8)  are  classical 
linear  terms,  while  the  final  two  are  nonlinear. 
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where  we  have  defined: 
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The  lift  has  been  adjusted  for  the  "leading  edge  force” 
in  the  classical  manner,  and  the  moment  is  about  the 
semi-chord.  If  (10)  and  ill)  are  expanded,  a  few 
important  points  become  evident: 


1 )  Linear  terms  are  recovered. 

2)  Between  the  prescribed  modes  of  motion  h(t)  and 
aft)  and  their  time  derivatives,  cross-coupling 
becomes  very  important. 

3)  No  single  term  requires  more  than  a  double 
quadrature  over  the  surface  of  the  airfoil,  and 
most  integrations  require  only  a  single  sweep. 


Reference  8  gives  a  full  explication  of  the  terms. 


Wagner’s  Nonlinear  Integral  Equation 

Wagner’s  integral  equation. shown  in  equation 
(13).  may  be  described  as  the  relationship  between  the 
quasi-static  circulation  T,,  (t)  bound  to  the  airfoil  and  the 
strength  of  the  trailing  vorticity  at  a  position  £  in  the 
wake.  The  extension  to  the  nonlinear  regime  effectively 
frees  the  wake  to  convect  under  local  velocities  by 
allowing  the  streaming  variable  s  to  integrate  over 
complex  space  where  the  wake  now  resides.  The 
resulting  equation  ( 14)  is  valid  when  we  are  careful  not 
to  integrate  from  one  branch  to  another,  as  proven  in 
appendix  A  of  reference  8. 


Typical  results  for  the  program  are  given  in  Figures 
4  through  7.  In  Figure  4  we  recreate  Karman’s  starting 
problem  by  imposing  a  step  function  of  magnitude  one 
degree  in  angle  of  attack  at  t=5  in  order  to  apply 
instantaneously  a  constant  bound  circulation  strength 
to  the  airfoil.  Figure  4a  is  a  trace  of  the  input  function, 
and  Figure  4b  is  a  picture  of  the  airfoil  and  wake  taken 
at  t=80.  The  free  stream  velocity,  which  is  the  same  for 
all  cases  presented  in  this  section,  is  equal  to  c/4  per 
time  step.  Figure  4b  shows  the  characteristic  roll-up  of 
the  starting  vortex  as  well  as  the  bowing  and  overall 
downward  displacement  of  the  wake. 


In  equation  (14),  the  wake  vortex  location  is 

measured  in  the  coordinate  system  fixed  to  the  plate, 
and  xmax  't>  is  taken  to  be  the  end  of  the  wake  farthest 
from  the  trailing  edge  xt  e  .  Equation  (14)  is  written  to 
take  advantage  of  the  discrete  vortex  method  of 
modeling  the  wake,  and  the  integration  relies  on  the 
fact  that  the  flow  is  assumed  to  be  non-dissipative.  This 
assumption  assures  that  an  effective  "history"  is 
maintained  in  the  wake.  It  also  guarantees  that  there 
are  two  unknowns  at  each  time  step:  the  strength  and 
the  position  of  the  vortex  formed  closest  to  the  trailing 
edge.  As  is  also  discussed  by  Scott,8  an  effective  and 
useful  solution  to  this  problem  is  to  use  an  offset  or 
"threshold"  distance  that  is  small  in  relation  to  the 
chord  length  (typically  c  .12  to  c  4).  Such  a  preset 
position  yields  Yu;  (s.ti  as  the  only  unknown.  P,,  if)  is  a 
known  quantity  at  each  time  step  because  it  depends 
only  on  the  velocity  and  position  of  the  airfoil,  which  in 
the  present  case  is  a  prescribed  quantity. 


The  preceding  analysis  has  been  incorporated  into  a 
fast  and  efficient  interactive  computer  program.  Inputs 
to  the  program  include  the  threshold  vortex  offset  dis¬ 
tance  and  the  history  of  the  motion  of  the  airfoil. 
Analytic  or  continuous  functions  are  not  required  as 
inputs,  so  the  user  may  create  effective  histories  for  each 
of  the  independent  and  overlapping  control  variables 
h’tt  and  cut/  by  construction,  using  pieces  of  functions 
such  as  ramps,  sinusoids,  and  steps  as  building  blocks. 
Our  experience  shows  that  the  only  limiting  parameter 
to  this  representation  is  the  commonly  seen  breakdown 
of  a  wake  formed  of  point  vortices.  Many  models  such  as 
the  Rankine  vortex,  the  cloud-in-cell  formulation.  ( >  and 
the  scheme  proposed  by  Mook  et  al.'2  have  been  offered 
as  viable  solutions  to  this  inherent  instability  and  will 
not  be  discussed  here.  Present  coding  allows 
considerable  leeway  in  the  choice  of  time  increment  and 
reduced  frequency  to  obtain  acceptable  results. 4 


In  Figure  4c,  the  lift  coefficient  asymptotes  to  its 
steady  value  as  the  starting  vortex  convects  down¬ 
stream.  In  addition,  Karman’s  starting  value  of  the  lift 
coefficient  is  recreated.  Just  as  in  the  classical  case, 
Ci(t)  attains  half  its  steady  magnitude  immediately 
after  the  motion  is  completed.  The  nose-down  pitching 
moment  (Figure  4d)  about  the  midchord  shows  initial 
overshoot  of  the  steady  value,  but  like  the  lift 
coefficient,  the  moment  coefficient  eventually  decays 
toward  the  steady  asymptotic  solution.  The  overshoot 
gives  insight  for  the  case  where  the  airfoil  is  uncon¬ 
strained  in  motion,  since  the  increased  nose-up  moment 
would  tend  to  increase  the  angle  of  attack  of  an  airfoil 
responding  to  induced  forces. 

Figure  5  shows  the  results  of  a  large-amplitude 
oscillation  of  the  height  of  the  airfoil  around  the 
reference  axis.  The  plunging  airfoil  experiences  the 
motion  shown  in  Figure  5a,  and  the  resulting  wake 
shows  the  large-scale  roll-up  patterns  of  Figure  5b.  This 
wavy  wake  resembles  an  early  stage  in  the  development 
of  the  Kdrmdn  vortex  street  as  exemplified  in  Van 
Dykel3  or  Mook  et  al.12  Figures  5c  and  5d  detail  the 
lead-lag  characteristics  of  the  motion,  with  Figure  5c 
showing  that  for  a  reduced  frequency  of  0.05  cycle  based 
on  the  semi-chord,  the  unsteady  total  lift  leads  the 
quasi-static  lift  by  approximately  29  degrees.  The 
amount  of  lead  or  lag  that  is  presently  predicted  by  the 
program  for  a  range  of  reduced  frequencies  from  0.005 
cycle  to  0.1  cycle  shows  divergence  from  the  classical 
results  for  linearized  small-amplitude  unsteady  motion 
derived  by  von  K&rm&n  and  Sears.3  However,  the 
amount  of  phase  divergence  is  fairly  small,  ranging 
from  six  degrees  at  low  reduced  frequencies  to  30 
degrees  at  the  higher  end,  and  is  purely  a  numerical 
error.  The  present  theory  may  be  rewritten,  as  outlined 
by  Scott8  and  McCuneS.  so  that  this  computational 
inaccuracy  in  calculating  the  normal  force  and  moment 
on  the  airfoil  is  eliminated. 

Figures  6  and  7  show  two  simple  maneuvers  and 
their  effects  on  the  loading  of  the  airfoil.  In  the  first 
maneuver,  the  airfoil  undergoes  a  sinusoidal  oscillation 
in  angle  of  attack  for  one  complete  cycle,  then  stops  its 
motion  altogether,  as  shown  in  Figure  6a.  Figure  6b 
shows  the  wake  curling  up  downstream,  while  Figure  6c 
shows  the  effect  of  the  maneuver  on  the  total  lift  coeffi¬ 
cient.  The  lift  vector  is  shifted  away  from  the  forcing 
function  a (t).  In  addition,  both  the  lift  and  the  moment 
(Figure  6d)  show  the  residual  effect  of  the  wake’s 
presence  after  the  maneuver  has  stopped.  The  overshoot 


in  both  cases  is  only  5  to  10  percent  of  the  maximum 
values,  and  the  sense  of  direction  of  the  overshoot  is  that 
of  the  last  completed  maneuver.  In  other  words,  the 
motion  of  the  airfoil  is  an  upward  pitch  as  it  completes 
the  sinusoid;  this  motion  sheds  counterclockwise 
vorticity  into  the  wake,  in  turn  creating  a  downwash  on 
the  airfoil.  The  downwash  causes  a  decrease  in  the 
airfoil's  lift  as  well  as  a  sizable  nose-up  moment. 

The  controlling  factor  for  the  amount  of  overshoot 
appears  to  be  the  magnitude  of  the  derivative  of  the 
forcing  function.  This  is  illustrated  in  Figures  7,  where 
the  input  is  a  step  function  in  angle  of  attack  followed  by 
the  negative  of  this  step  30  time  increments  later 
(Figure  7a).  Figure  7b  shows  that  the  wake  has  a 
"square  wave”  appearance,  and  the  lift  augmentation 
(Figure  7c)  is  instantaneously  on  the  order  of  20  percent 
of  the  quasi-static  value.  Step  changes  cause  the  airfoil 
to  move  at  high  velocities  (computationally  finite  but 
theoretically  infinite).  These  velocities  then  introduce 
large  amounts  of  vorticity  into  the  wake.  In  the  present 
study,  instantaneous  lift  augmentation  on  the  order  of 
50  percent  has  been  seen.  And,  as  shown  in  Figure  7d, 
the  moment  augmentation  for  the  test  case  is  sub¬ 


stantially  larger  than  that  observed  in  the  preceding  Cl 
history. 

Conclusion 


The  strong  points  of  the  present  computational 
method  are  its  speed  and  flexibility.  Since  we  solve  for 
only  one  unknown  per  time  step,  the  results  may  be 
displayed  interactively  on  a  graphics  screen.  A  typical 
run  of  80  time  steps,  which  may  extend  the  wake  20 
chord  lengths  or  more  downstream,  takes  only  5  to  10 
minutes  to  generate  when  run  on  a  DEC  VAX  11-750 
computer.  In  fact,  the  speed  of  the  graphics  plotting 
package  is  usually  the  limiting  factor  in  such  a 
calculation.  In  terms  of  flexibility,  the  program  can 
accept  any  combination  of  sinusoids,  pieces  of  sinusoids, 
step  functions,  or  constant  values  in  each  of  the  input 
variables  h(t)  and  aifA  The  user  is  free  to  "fly”  a  given 
maneuver  on  his  own,  having  only  to  "construct"  a 
realistic  forcing  function  history  with  an  interactive 
input.  The  advantages  of  the  method  are  therefore 
threefold:  the  user  input  is  flexible,  the  program  is  fast, 
and  the  results  are  completely  nonlinear  and  exact. 


Figure  4.  The  starting  problem:  Instantaneous  loading 
of  an  airfoil. 
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Figure  4a.  Input  angle  of  attack  history  taith. 
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Figure  4b.  Airfoil  and  wake  at  t=30. 
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Figure  5.  Sinusoidal  oscillation:  Time  variant  loading 
of  an  airfoil. 
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Figure  5a.  Input  plunging  velocity  history  (h lt)i. 
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Figure  5b.  Airfoil  and  wake  at  t=80. 
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Figure  5c.  History  of  lift  coefficient. 

Figure  4c.  History  of  lift  coefficient. 
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Figures.  Sinusoidal  maneuver:  Time  variant  loading  Figure  7.  Completed  step  function:  Instantaneous 

and  unloading  of  an  airfoil.  loading  and  unloading  of  an  airfoil. 


Figure  6a.  Input  angle  of  attack  history  'a utJi 


Figure  7a.  Input  angle  of  attack  history  'o tit)). 


Figure  6c.  History  of  lift  coefficient. 


Figure  7c.  History  of  lift  coefficient. 


Figure  6d.  History  of  moment  coefficient. 


Figure  7d.  History  of  moment  coefficient. 
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INTERACTIVE  AERODYNAMICS  OF  WINGS  IN  SEVERE 

MANEUVER*1 

J.E.  McCUNE 

Massachusetts  Institute  of  Technology,  Cambridge,  Massachusetts 


ABSTRACT 

Thi«  paper  describes  an  interactive  technique  for  the 
study  and  analysis  of  the  aerodynamic*  of  wing*  in  *evere 
maneuver.  Both  the  2D  airfoil  case  and  the  cross-flow  prop¬ 
erties  for  slender  3D  wings  are  addressed.  To  begin  the 
study,  the  non-linear  interaction  of  a  2D  airfoil  undergo¬ 
ing  large-amplitude  unsteady  motion  with  its  wake  of  shed 
vorticity  was  studied  interactively  on  the  computer  using 
quick  and  efficient  codes.  The  user  can  now  input  a  ‘ma¬ 
neuver”  and  study  and  observe  on-line  the  non-linear  wake 
evolution  and  airfoil  response.  The  allowable  maneuvers 
for  the  airfoil  presently  include  any  combination  of  pitch¬ 
ing  and  plunging,  suddenly  imposed.  Also,  airfoil  response 
to  sudden  large- amplitude  gusts  of  any  shape  and  relative 
passage  speed  can  be  observed  and  analyzed.  A  detailed 
review  of  progress  and  results  from  this  part  of  our  work 
is  provided  by  Scott. [l]  In  the  present  report  we  empha¬ 
size  application  of  related  techniques  to  3D  slender  wings. 
We  also  provide  information  on  improved  methods  of  cal¬ 
culating  lift  and  moment  via  wake  integrals,  emphasizing 
conservation  of  impulse,  for  the  2D  as  well  as  the  3D  case. 

1  Introduction 


role  of  this  vorticity,  emanating  via  "leading-edge  separa¬ 
tion,”  is  to  maintain  acceptably  smooth  flow  at  the  wing 
edges.  At  sufficiently  large  Reynolds  numbers,  this  smooth 
flow  requirement  is  often  called,  even  in  this  context,  a 
“Kutta  condition.” 

The  present  paper  emphasises  recent  progress  in  the  de¬ 
velopment  of  a  new  method  for  studying  on-line  the  dynam¬ 
ical  effects  of  these  two  cross-flow  wakes  for  the  3D  wing. 
The  method  proceeds  in  analogy  with,  but  also  extends, 
techniques  used  in  the  2D  airfoil  case.  Additional  ideas 
and  theoretical  framework  needed  for  the  3D  case  are  out¬ 
lined.  It  is  shown  that  the  wake  structure*  above  a  delta 
can  be  determined  in  a  manner  similar  to  the  2D  airfoil 
problem,  and  that  these  structures  have  a  history  unique 
to  any  given  maneuver. 

As  discussed  in  the  body  of  the  paper,  the  model  of  the 
wake  structure  used  is  somewhat  idealized,  in  that  the  de¬ 
velopment  during  wake  roll-up  of  vortex  “cores”  at  finite 
Reynolds  Numbers,  with  attendant  regions  of  distributed 
total  pressure  loss,  is  not  addressed.  Similarly,  the  related 
occurrence  of  vortex  breakup  is  not  included  in  the  present 
model.  Our  recent  efforts  to  deal  with  these  matters,  as 
an  elaboration  of  the  technique  described  here,  will  be  re¬ 
ported  elsewhere. 


One  important  use  of  the  non-linear  2D  airfoil  study  is 
to  provide  a  test  of  the  exact  wake  evolution  method,  ap¬ 
plied  at  low  amplitudes,  against  classical  linearized  airfoil 
theory.  This  tests  the  accuracy  and  efficiency  of  the  present 
approach,  at  least  in  the  classical  limit.  But,  in  addition, 
the  2D  work  has  set  the  stage  for  significant  advances  in  the 
study  of  the  3D  non-linear  aerodynamic  response  of  low-to- 
moderate  A  wings  (delta  and  other)  to  severe  imposed 
maneuver*. 


2  Non-Linear  Unsteady  Wake  Convec¬ 
tion 

2.1  2D  Case 

Any  free  vorticity  associated  with  unsteady  plane 
2D  fluid  motion  in  the  incompressible  limit  obeys  the 
Helmholtz  relation 


It  has  long  been  recognized  that  at  Urge  but  finite 
Reynolds  Nr*,  the  cross-flow  at  any  chordwise  station  of 
a  slender  wing  must  include  two  “wakes*  (actually,  wake 
traces)  representing  concentrated  vorticity,  often  partially 
relied  into  “cores,”  convecting  above  the  wing  surface.  The 

'This  work  is  dedicated  to  my  teacher.  Prof.  Wm.  R.  Sears, 
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ported  by  NASA  Lanffley  under  Grant  Nr.  NAC-i-65*. 
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where  V  is  the  planar  fluid  veloctiy,  V  =  (u,t>,0),  and 
0  =  curl  V  =  (0, 0,  n,(x,  y,  t)) .  Formally,  therefore,  if 
A(1)  =  A and  A{1)  =  n(z,y, /)  are  two  independent 
characteristic  solutions  of 


(2.2) 


1 


with  appropriate  boundary  (/  initial)  condition*,  then  fl«  = 
fi,(n,A).  Boundary  condition*  can  b*  *p«cified  *uch  that 
one  of  the**  characterise  variable*,  “n,*  say,  i*  constant 
on  streakline*  drawn  at  fixed  t.  The  remaining  variable,  A, 
can  be  choeen  to  be 

A  =  t  —  r  (2.3) 

where 


Thu*,  in  the  20  caee  we  can  write  A£  =  A$(A),  and  d r  = 
dT(A)  a  -  =  P(A)dA  where 

A  5f'?  (2.9) 

and 


—  =  _  +  <  V  >  V7=  1  (2.10) 

Dt  dt 


(2.4) 


“r"  is  then  the  Eulerian  “drift  time*  with  the  usual  sim¬ 
ple  physical  interpretation  [2].  But  if  the  velocity  field  i* 
discontinuous  a*  in  the  classical  model  of  concentrated  vor¬ 
tex  sheet*  representing  the  wakes  behind  wings  and  airfoils, 
this  formulation  may  not  be  convenient. 

It  turn*  out,  however,  that  when  a  concentrated  vortex 
sheet  is  embedded  in  the  flow  the  variable  v  can  be  con¬ 
veniently  replaced  by  a  continuous  function,  say,  which 
instead  of  (2.2)  satisfies 


=  ^+<V  >.V„e=0 
AV  Vne  =  0 


i.e.,  ije  i»  constant  on  the  sheet,  and  the  sheet  is  a  streak 
surface  (streakline  in  2D),  whose  deformation  is  controlled 
by  convection  at  the  mean  velocity,  <  £  >: 


<  £  >3  ^  (K+  ^  K')  (2-6) 

As  with  a  vortex  sheet, 


AV=V+-r  (2.7) 

r.as  no  component  perpendicular  to  the  sheet  at  any  time, 
as  stated  in  (2.3(b)). 


Here,  7  is  the  “mean*  drift  time  of  fluid  elements  within 
the  vortex  sheet. 

The  variable  rje  can  be  expressed  in  the  form  Vc  — 
y  -  y<(7,  t)  where  =  v.  With  r»e=const.  on  the  streak¬ 
line,  we  have  y  =  const. +  Yc(7,t)  on  that  line.  When  7  is 
expressed  in  term*  of  x  at  any  t  we  define  a  space  curve, 
y  =  y(r;t),  giving  the  usual  instantaneous  picture  of  an 
airfoil  wake.  (See  Figures  1  and  2.) 

A*  noted  above,  condition  (2.8)  was  derived  by  requir¬ 
ing  that  there  be  no  net  force  anywhere  on  the  free  vortex 
sheet.  For  systems  of  conserved  global  circulation,  this  cor¬ 
responds  to  guaranteeing  conservation  of  impulse  for  the 
airfoil-plu*-wake. 

In  the  linear  case,  in  which  the  wake  is  assumed  simply 
to  lie  along  the  2-axis  extending  from  the  trailing  edge  with 
the  vorticity  within  it  convected  at  free-stream  speed  Uoo, 
the  above  exact  results  reduce  to 


dA  1 

d  x  Uoo 

so  that,  if  -7(2,  t)  is  the  wake  vorticity  in  that  limit, 

dr(X)(l‘near)  _  ^xJ-jdx 

and 


In  a  distance  da  along  the  sheet  at  any  fixed  t  the  ele¬ 
mental  change  in  circulation,  dr,  is  determined  by  the  local 

jump  across  the  sheet  of  the  velocity  potential  <p{x,  y,  f)  such  dP  in^  dA  1  d  /r-(lin)  /,  *  -  Zo  \  'l  _  /,  z~ 

: hat  i\  dx^  UvodtV  V  U„  ) )  7V  V, 


-dr  =  d$A~~  =  d»~ Ad  =  dA<fi  (fixed  t) 
da  da 

In  the  2D  unsteady  case  the  strength  of  this  potential  jump, 
A <p,  occurring  at  the  sheet  specified  by  =  constant,  is  a 
function  of  both  x  and  t  which  must  be  determined.  But 
application  of  the  Bernoulli  equation  on  either  side  of  the 
sheet,  together  with  the  requirement  that  there  be  no  jump 
in  static  pressure  across  the  sheet,  yields  the  restriction 

^+<V>VA*  =  0.  (2.8) 


reproducing  the  classical  result  in  the  linear  limit. 

In  the  actual  non-linear  situation  A  replaces  z  as  the  use¬ 
ful  parameter  along  the  vortex  sheet  representing  the  wake, 
whether  the  wake  is  distorting  and  rolling-up  or  not.  Thus, 
for  example,  the  total  wake  circulation  at  any  given  time  t 

is 

rwake(0=  /  rfAr'(I).  (2.11) 

j  0 

Further,  the  (plane,  2D)  velocity  induced  by  the  wake  can 
be  written  in  complex-variable  notation  a* 


uwaks(*’  V'  0  —  ,uwake(*’  1'  *)  —  wvt*kt(*'  *) 


(M2) 


Here,  /  =  *  -My,  and  *,  =  *.  -My,  is  the  complex  location 
at  t  of  the  vortex  element  having  circulation  dr(A). 

Thie  formulation  provides  a  powerful  tool  for  writing 
down  exactly  all  the  various  non-linear  effects  related  to 
the  wake  in  a  manner  which  conveniently  follows  the  con¬ 
tortions  of  the  wake,  including  roll-up,  as  it  evolves.  Since 
A  is  conserved  during  the  history  of  each  fluid  element  in 
the  wake,  it  provides  a  convenient  label  for  each  such  el¬ 
ement.  In  the  present  example  A  is  the  (earlier)  time  at 
which  a  given  element  entered  the  wake  at  the  airfoil  trail¬ 
ing  edge.  Figure  1  illustrates  a  portion  of  the  evolution  of  a 
streakline  constituting  a  vortex  sheet,  and  Figure  2  shows 
a  typical  example  of  a  computer-generated  version  of  the 
wake  behind  an  oscillating  airfoil. 

On  the  computer,  of  course,  the  above  continuous  vortex 
sheet  model  is  replaced  by  an  approximate  discretized  ver¬ 
sion.  Individual  discreet  vortices  of  strength  6T(A)  replace 
the  line  elements  of  the  sheet.  Details  of  how  to  do  this 
properly  so  as  to  maximize  the  effectiveness  of  the  approx¬ 
imation  have  been  discussed  by  Mook,  et.  al.  [3|  and  more 
recently  by  Scott  [lj. 

The  computer  screen  provides  an  almost  ideal  tool  for 
depicting  and  understanding  wake  evolution  in  this  frame¬ 
work.  At  each  time  step  each  fluid  element  (and  free  circu¬ 
lation  element)  is  advanced  to  its  new  location  in  z  and  y, 
moving  in  the  field  of  all  its  neighbors.  Since  A  is  constant 
for  any  fluid  element  and  t  is  known,  then  7  is  determined 
at  each  new  location  at  each  t.  In  effect,  the  computer  au¬ 
tomatically  maps  the  drift  time  field.  When  all  the  fluid 
elements  in  a  given  wake  are  shown  in  their  new  positions 
at  each  new  time,  t,  for  example,  the  line  through  them 
and  the  trailing  edge  at  that  moment  is  a  streakline  which 
depicts  the  momentary  shape  and  location  of  the  vortical 
■wake.  We  discuss  in  Section  3  how  the  individual  strengths, 
dT,  of  the  circulation  elements  are  determined  (the  Wagner 
equation). 


Although  this  model  is  certainly  oversimplified  because 
of  such  effects  as  separation  of  the  return  boundary  layers 
on  the  wing  upper  surface  (leading  to  ‘secondary*  and  ‘ter¬ 
tiary*  vortices  [9|,  (10|),  it  offers  an  important  improvement 
for  real  flows  over  the  completely  inviscid  claesical  theory 
(Re  5  oo)  with  its  singular  flow  at  the  edges  [11).  This 
postulated  wake  pattern  is  often  referred  to  as  the  result 
of  ‘leading-edge  separation*  in  the  present  context.  The 
strength  of  the  vorticity  is  presumed  to  be  determined  so 
as  to  provide  Kutta  conditions  (smooth  flow)  at  both  wing 
edges.  (Section  3.) 

At  higher  angles  of  attack,  even  in  steady  flow,  the  model 
pattern  in  Figure  3  has  a  further  defect.  Vortex  breakup 
[12] ,  [13],  [U|  can  be  expected  to  occur  at  some  conditions 
of  interest  with  a  consequent  sudden  and  intense  change 
at  some  chordwise  location  of  the  vorticity  pattern  above 
the  wing.  This  paper  does  not  attempt  to  deal  with  vortex 
breakup  phenomena. 

In  the  unsteady  case,  especially  for  3D  wings  in  severe 
maneuver,  the  classical  wake  pattern  also  has  other  de¬ 
fects.  Clearly,  for  example,  a  pilot  could  well  be  able  to 
roll  or  pitch  so  rapidly  as  to  catch  and  move  through  his 
own  wake  pattern.  Understanding  wing  response  to  that 
sort  of  condition  will  require  improved  understanding  of  the 
vorticity  distribution  within  the  cores  themselves  and  any 
associated  total  head  loss  within  the  wake  system.  [IS],  [16] 
For  that  matter,  the  same  improved  understanding  of  the 
cores  and  total  pressure  defects  appears  to  be  required  to 
interpret  vortex  breakup  jl2].  In  the  following  discussions 
we  largely  ignore  such  defects,  leaving  their  treatment  for 
future  elaborations  of  the  technique  we  wish  to  describe 
here. 

The  evolution  of  the  free  3D  wake  pattern  over  the  wing 
proceeds  in  a  manner  analogous  to  the  2D  case,  complicated 
only  by  the  need  for  an  additional  characteristic  variable 
in  the  chordwise  direction.  To  maintain  zero  net  force  on 
the  free  wake  system  each  circulation  element,  -  Ad(z,  z,  f ), 
once  it  enters  the  wake  pattern  at  a  wing  edge  location,  is 
convected  unchanged  at  the  speed  determined  by  the  mean 
velocity  at  the  sheet,  including  the  motion  “induced”  by  all 
its  neighboring  vortex  elements.  Thus,  again, 


2.2  3D  Case 


<  V  >  =  0  (2.1*3) 


The  classical  picture  of  the  3D  vortex-wake  pattern  above 
a  delta  wing  in  steady  flow  at  large  Reynolds  Nr.  and  mod¬ 
est  angle  of  attack  is  shown  in  Figures  3  and  4.  At  any 
chordwise  station  the  trace  of  the  double  wake  (somewhat 
idealized)  is  as  illustrated  in  the  inset  of  Figure  3.  The 
necessity  for  a  wake  pattern  resembling  this  structure,  in 
order  to  provide  acceptably  smooth  flow  at  the  wing  edges 
at  large  Re,  has  been  recognised  by  many  authors  [4],  [S], 
[Sj,  [7],  [8|.  Even  though  the  vorticity  may  roll  up  partially 
into  two  strong  ‘cores,*  the  strength  of  these  cores  must 
vary  in  the  chordwise  direction,  which  requires  the  exis¬ 
tence  of  the  joining  umbilical  vortical  sheets  from  the  edges 
to  the  cores,  as  illustrated. 


where  V  =  (u,u,u>)  and 

<£>=j(*r+)n  (2.H) 

Once  more  we  define  the  mean  drift  time,  7,  now  in  the 
3D  field,  such  that 

^+<V>  V?=l  (2.15) 

and  reintroduce  the  characteristic  variable  (now  3D) 

A  =  t  -  f  (2.16) 
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which  i*  constant,  u  before,  for  each  fluid  element  moving 
in  the  vortex  wake  eyetem- 

What  emerge*  i*  the  deecription  of  the  overall  3D  wake 
pattern,  however  contorted,  made  up  of  twitted  etreaklinet 
of  vorticity  extending  at  any  inetant  from  the  plane  of  obeer- 
vation  at  a  given  chordwiee  location  to  the  upetream  wing 
edge  location  at  which  the  corresponding  vortex  elemente 
have  entered  or  are  entering  the  wake  tyitem.  To  detcribe 
the  shape  and  location  of  these  streaklines  we  again  need 
the  characteristic,  ue,  as  defined  in  (2.5),  and  one  additional 
characteristic 

f.  =  z  -  Ze{z,y,z,t)  (2.17) 

where  _ 

f  5|+<v>^  =  ,  (2.18) 

For  the  latter  we  can  choose  the  boundary  condition  Zc  =  0 
at  or  along  any  wing  edge  so  that  Ze  is  just  the  chordwiee 
distance  at  any  t  and  x  a  given  element  of  the  wake  has  trav¬ 
elled  since  entering  the  system  at  the  upstream  chordwiee 
location  and  at  a  corresponding  earlier  time  A.  Thus, 

=  0  for  any  wake  element  entering  the  pattern  at  the 
wing  vertex,  and  U-*  for  any  element  just  entering  at  the 
plane  of  observation  at  x.  The  tj»  and  ^  variables  define 
instantaneous  surfaces  in  the  3D  case,  and  their  intersec¬ 
tions  are  the  streaklines  making  up  the  wake  pattern  at  any 
given  moment. 

In  analogy  with  the  2D  case  one  can  express  the  variables 
vc  and  u  in  the  forms 

ic  =  y  -  Ye(f,S',t) 

U  =  *  -  Zc(7,1c,t) 


f,  in  turn,  is  a  function  at  each  z  and  t  of  z,  the  height  of 
the  wake  element  above  the  wing  surface.  For  example,  in 
the  linear  model  of  a  wake,  with  wake  elements  conveeted 
:n  straight  lines  at  speed  C/oo  a*  depicted  in  Figure  4(b), 
we  have 

?  =  UZ  sina  1  f*  =  Zed««  1  ”e  =  Ked«e 

y  =  ^edge(7edge)  ^  /U^fcota 
z  =  Zedge  +  UaoCO»aJ 


so  that  the  simple  surfaces 

2  “  2edge 
y  ~  ^edge 


h 

tan  a 
[h 
tan  a 


plane.  In  the  unsteady  case,  these  traces  will  deform,  con¬ 
tort,  and  appear  to  roll-up,  much  as  if  we  were  observing  a 
2D  airfoil  with  two  wakes  rather  than  one. 

In  Figure  5  we  illustrate  the  sort  of  wake  surfaces,  and 
their  traces,  to  be  expected  above  a  delta  wing  undergoing 
yaw  and  roll.  Simultaneous  observation  of  the  wake  trace 
behavior  at  several  crose-flow  planes  can  provide  significant 
insight  as  to  the  interactive  aerodynamics  associated  with 
violent  maneuvers.  On-line  use  of  computer  graphics,  in¬ 
cluding  split  screen  capabilities,  promises  to  be  very  helpful 
in  understanding  the  large-scale  unsteady  behavior  of  the 
wing  and  its  wake. 

On  taking  account  of  the  right  hand  rule  in  the  definition 
of  circulation,  each  circulation  element  in  each  trace  at  fixed 
z  and  t  can  now  be  described  in  terms  of  the  above  variables 

as 


tfroxr  _  +dAi(i',\)roKr 


in  view  of  (2.13).  (In  the  steady  case,  with  nothing  changing 
in  time  at  a  given  z,  only  the  first  variable  need  appear  and 
dr  =  dr(*  -  Zc(z,y,z))  =  r'(f.)df.)  in  each  wake. 

The  drift  time  f  observed  at  z  and  t  for  any  element  in 
either  of  the  two  wake  traces  is  related  to  fe,  so  that  we 
can  find  A(f,-,s,x,t),  and/or  the  inverse.  Then,  the  total 
circulation  at  fixed  z  and  t  in  either  trace  (starboard  or 
port)  can  be  written 


rwakeD'°*r(2-‘)  =  /  dr(f<J)(f< 
walte  J(  fixed*..) 


(STflD.POftr) 
/**  ^STBD.POHr 

■/, d<— — 

(dA4>  <?Ad  dA\ 
V.  v  3z  dt  dS'J 


STBD.rORT 


Note  that  in  the  linear  case  or  for  slender  wings  at  mean 


angle  of  attack  a, 


As  discussed  in  Section  3.2,  the  incremental  circulation 
associated  with  each  vortex  element  is  determined  just  as 
it  enters  the  wake  pattern  so  as  to  always  maintain  the  pos¬ 
tulated  smooth  flow  at  the  edges.  After  that,  each  element 
is  simply  conveeted,  at  fixed  elemental  circulation,  to  new 
locations  in  the  wake  pattern. 


3  Determining  the  Vortex  Strength 
Distribution  in  the  Wakes 


are  formed.  The  more  realistic  non-linear  case  follows  sim¬ 
ilar  rules  and  is  depicted  in  Figure  4(a). 

When  observed  at  a  given  chordwiee  location,  z,  in  the 
corresponding  cross-flow  plane  including  the  wing  trace,  the 
vortex  pattern  over  the  wing  create*  a  trace  (actually,  a  pair 
:f  traces)  of  the  wake  a*  it  passe*  through  the  observation 


3.1  2D  Airfoil  Case 

In  our  treatment  of  this  case  we  assume  the  Kutta  con¬ 
dition  to  hold  at  the  airfoil  trailing  edge  at  any  instant. 
The  airfoil  itself  is  assumed  to  be  a  fiat  plate1  in  arbitrary 

‘a  central  2D  airfoil  profile  can  be  treated  in  tbe  earae  man- 
ner.  with  only  *  minor  elaboration  of  the  technique  used  here 


«.  s 


plunging  and/or  pitching  motion. 

A*  already  discussed,  the  uneteady  airfoil  motion  requires 
the  existence  of  a  wake  such  as  illustrated  in  Figure  2,  so 
that  the  net  circulation  of  the  wake-plus-airfoil  is  constant. 
No  linearisations  of  the  airfoil  motion  of  wake  contortions 
are  admitted  in  the  present  treatment. 

The  airfoil  circulation  can  be  described  conveniently  in 
two  parts.  The  first,  r0(f).  is  the  “quasi-steady*  Kutta 
value  associated  with  the  airfoil  motion  itself  and  calculated 
as  if  no  wake  were  present.  Thus,  r0(t)  is  known  for  any 
specified  airfoil  motion.  The  general  expression  is 


To  =  -wcUao  sin  o(f)  +  2  /  diio{t,t)y 


The  quantity  dwake(2,f)  can  obtained  readily  from 
(2.12).  Using  a  system  i  =  e_,a(,,(*  -  »h(£))  =  2  +  ty  ro¬ 
tating  and  plunging  with  the  airfoil,  the  normal  component 
of  the  apparent  wake  velocity  at  the  airfoil  is  determined. 
Inserting  the  result  in  (3.7)  yields  an  expression  for  Ti(f) 
in  terms  of  an  integral  over  the  wake  vorticity.  Details  of 
this  calculation  are  available  in  Eef.  [1|. 

The  next  step  is  to  apply  the  Kelvin  Theorem  in  the  form 


where  *c*  is  the  airfoil  chord,  a(t)  is  the  angle  of  attack, 
i  is  the  instantaneous  chordwise  coordinate  along  the  air¬ 
foil,  and  t>o(2,t)  is  the  actual  normal  component  of  fluid 
velocity  at  the  oscillating  and  plunging  plate  necessary  to 
accommodate  its  unsteady  motion.  (See  Figure  6).  Then, 
without  approximation,  (see  Ref.  (l|  for  more  details) 

Oo  =  ficoia(t)  —  2d  (3.2) 


Jf  c- 

r0(f)  =  -ire(l/ooSina  -  A  cos  a) - —  a.  (3,3) 

The  bound  vorticity  associated  with  To  and  Oo  is  called 
7o(i,  £)  and 

r0(£)  =  f  70  (z,t)dx.  ( 3A ) 

J'i 


Defining  the  angle  variable  0  such  that  2  =  |cos/?,  we 


7o(2,  t)  =  -2  (C/«,sina  -  Acoea) 

-  ca  sin/? 


1  -  cos  0 


This  “quasi-steady*  bound  vorticity  would  be  all  that  is 
present  if  it  weren’t  for  the  wake.  But,  as  illustrated  in 
Figure  6,  and  calculated  in  Eq.  (2.12),  the  wake  causes 
an  apparent  upwash  at  the  airfoil  which  must  be  cancelled 
cut  by  the  action  of  additional  bound  vorticity  on  the  air¬ 
foil,  71  (2, f).  The  amount  needed  is  given  by  the  theory  of 
conjugate  functions  in  the  form 


„  1-costf /'*(-2*wlk#)(I  +  cosr) 

— cY.r-ro.5 —  (3-6) 


where,  as  in  (3.5),  the  Kutta  condition  has  been  applied  at 
the  trailing  edge.  Correspondingly,  the  additional  circula¬ 
tion  on  the  airfoil,  ri(t),  is  given  by 


rairfoil(0  +  rwake(‘)  =  con*tMt  =  r(°)  (3-8) 

where  T( 0)  is  any  existing  initial  or  steady-state  airfoil  cir¬ 
culation.  For  example,  if  the  airfoil  maneuver  begins  from 
steady  flight  at  angle  of  attack  oo, 

T(0)  =  -ecsin  a0. 

But  we  have 

^airfoil  =  + 

and  rwake(f)  is  written  out  in  Eq.  (2.11),  also  in  terms 
of  an  integral  over  the  wake  circulation  elements.  Cancel¬ 
lations  occur,  and  the  final  result  is,  in  complex-variable 
notation, 


ru(0  -  r(o)  =  Re  dx r'(i),  Mri 

Jo  V 

where  “Re”  means  real  part  implied  and 

i,  =  e'“  (*.  -xh(t)) 

with  i,  as  defined  below  (2.12). 

Equation  (3.9)  was  first  derived  by  Wagner  [17),  in 
the  linearised  limit,  and  used  by  many  authors  to  under¬ 
stand  airfoil  response  to  sudden  starts,  sudden  (but  small) 
changes  in  angle  of  attack,  and  flying  through  gusts.  [I8i, 
[19|  It  seems  remarkable  that  the  non-linear  version  (3.9) 
of  Wagner’s  equation  is  essentially  the  tame  as  his  origi¬ 
nal  except  that  2.(A,f)  is  complex,  reflecting  the  distortion 
and  roll-up  of  the  wake.  The  principal  mathematical  nicety 
in  the  non-linear  case  is  the  need  to  determine  the  correct 
branch  of  the  complex  square  root  occurring  in  the  integral. 
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The  left- hind  tide  of  (3.9)  ie  known  as  a  function  of  t  for 
any  given  impoeed  airfoil  motion  (m«  Eqn.  (3.3),  for  exam¬ 
ple),  to  Eqn.  (3.9)  ie  an  integral  aquation  for  the  strength*, 
d r(I),  of  the  circulation  elemente  in  the  wake. 

Operationally,  on  the  computer,  uee  of  Eqn.  (3.9)  to 
determine  the  wake  vorticity  ie  relatively  etraightforward, 
except  for  finding  the  initial  incremental  vortex  element 
at  the  very  etart  of  a  maneuver.  Difficultiee  ariee  there 
becauee  the  kernel  on  the  tK*  of  (3.9)  it  tingular.  An  ef¬ 
fective  method  for  treating  thit  crucial  first  moment  of  the 
“starting  problem"  it  detcribed  in  Ref.  (1).  After  the  pro¬ 
cedure  hat  been  properly  ttarted,  however,  it  runt  easily 
becauee  once  the  atrength  of  a  vortex  element  hat  been  de¬ 
termined  by  (3.9),  at  the  inttant  it  enter*  the  wake,  that 
strength,  a*  measured  by  its  incremental  circulation,  re¬ 
mains  unchanged  at  the  element  convectt  to  its  subsequent 
positions  in  the  wake. 

3.2  3-D  Case.  Smooth  Flow  at  Both 

Edges 

From  Kelvin’s  theorem  we  know  the  net  circulation  in 
each  crose-flow  plane  it  zero: 

rjreo  +  feosT  +  Ftracb  —  0  (3.10) 

We  treat  the  3D  case  in  this  paper  in  the  slender-wing 
(low  /R)  limit,  determining  the  velocity  and  potential  field 
in  each  cross-flow  plane  on  a  quasi-two-dimensional  basis 
in  the  spirit  of  Ref.  [11].  In  that  limit  we  must  determine 
the  2D  potential  solution  in  each  observation  plane  corre¬ 
sponding  to  the  concentrated  vorticity  on  the  wing  trace 
and  in  the  traces  of  the  two  wakes  (Figure  7).3 

In  order  to  treat  a  general  large-amplitude  maneuver  and 
still  maintain  an  inertial  frame  "almoet"  moving  with  the 
wing,  we  specify  a  system  moving  with  the  mean  or  initial 
steady  flight  condition,  as  sketched  in  the  figure.  Thit  ini¬ 
tial  or  mean  forward  speed  is  !/<»,  with  angle  of  attack  5- 
A  cross-flow  observation  plane  at  any  fixed  z  (analogous  to 
the  computer  screen)  is  oriented  normal  to  the  wing  sur¬ 
face  in  its  initial  or  mean  configuration.  The  observer  in 
this  plane  tee*  an  incoming  steady  free  stream  of  speed  Uao 
at  angle  5  to  the  z  axis.  Additional  pitch,  yaw,  or  roll  of 
the  wing,  as  well  as  plunging  or  sideslip  of  the  wing  ver¬ 
tex,  are  perceived  in  each  observation  plane  as  vertical  or 
sideways  displacements,  Xo(z,t)  and  Ya(z,t).  Acceleration 
or  deceleration  of  the  vertex  in  the  z-direction  is  seen  as  a 

JF or  3D  wises  sod  their  wakes  durinf  sudden  sod  violent  ma¬ 
neuver  the  -slenderness  direction*  becomes  somewhat  ambifu- 
out-  Near  mod  on  the  low  wing,  (or  example,  quantities 
may  vary  slowly  io  spacs  only  in  ths  instantaneous  cbordwise 
direction,  whereas  this  is  eeosrally  misaliensd  with  ths  orieina! 
t-axis  of  Fieure  7.  Ths  maio  wahs  structure,  on  the  other  band, 
may  still  tend  to  be  ‘slender*  with  respect  to  the  orieina!  (licht 
alienments.  Ultimately,  of  course,  use  of  the  full  3D  relations 
is  requirsd.  Ws  ars  suec**tine  hers,  however,  that  this  be  done 
after  fully  exploitinc  the  additional  iasiebta  to  be  eained  in  the 
Jones-Munk  tradition. 


growth  or  shrinking  of  ths  wing  trace  at  any  s.  For  ma¬ 
neuver*  with  v*ry  rapid  forward  acceleration  or  complete 
changes  in  flight  direction,  a  frame  fixed  in  the  atmosphere 
may  be  preferable,  with  the  wing  flying  "through"  the  ob¬ 
servation  plans  or  planes.  Of  count,  one  could  also  choose 
an  observation  system  fixed  in  ths  wing,  with  attendant 
complication!  associated  with  non-inertial  frames  of  refer¬ 
ence.  For  our  present  purposes,  we  adopt  the  scheme  shown 
in  Figure  7. 

In  complex-variable  notation  the  quaai-2D  cross-flow  ve¬ 
locity  field  associated  with  the  wakes  ie,  with  Z  as  X  f  tF, 

^wake  ~  *^wake  =  ^wake(^) 

1 

2iri 

+ 

in  analogy  with  (2.12).  Here,  w*  have  used  (2.19)  and 
£/r*D(f«.f.»)  Z:°*T{u,t,z)  are  the  instantaneous 

complex  locations  in  any  given  croes-flow  plane  (x,  y)  of  the 
vorticity  element*  in  the  starboard  and  port  wake  traces  re¬ 
spectively.  In  the  following  we  also  have  uee  for 

^wakeW  —  t  ^wake 

“  ^wake  ~  ‘^wake • 

The  boundary  condition  at  the  wing  trace  a*  viewed  in 
each  cross- flow  plane  can  be  written  in  a  form  somewhat 
similar  to  that  in  (3.2)  for  the  2D  airfoil.  Using  the  coor¬ 
dinates  suggested  in  Figure  7  the  actual  component  of  ve¬ 
locity  in  the  (x,y)  plane  normal  to  the  instantaneous  wing 
trace  can  be  written 


f  r.  dr»r‘°(f.,A)  1 

l  Jo  U  du  Z-Zr,D 

f.  dr°»r(f., a)  i  \ 
J.  au  df.  z  -  zr**  { 


00{y,z,t) 


D,Xo(z,t)  is  ,  D,Yo(z,t) 
Dt - C0*K+— DT- 


a'in  K 


.  D,K{z,  t) 
y  Dt 


(3.12) 


where  For  a  low  /4i  wing  in  steady  forward 

Sight  without  pitch  or  yaw,  for  example,  W  =  t/oocosa, 
since  the  additional  fluid  speed  in  the  chordwise  direction 
is  negiibls  because  of  the  wing  slenderness.  During  pitch 
and  yaw,  however,  W  requires  adjustment  to  allow  for  the 
instantaneous  flight  angles  and  the  fact  that  the  chordwise 
(or  "slendernete")  direction  is  then  misaligned  with  the  z- 
axis  chosen  in  Figure  7. 

Starting  from  the  mean  or  initial  steady  Sight  condition 
described  above,  an  observer  on  the  wing  see*  an  effective 
crote-fiow  from  below  of  Woo  sin  5,  the  normal  component 
of  which  is  Uao  sin  5eos  K.  Eqn.  (3.12)  then  allow*  for  any 
additional  unsteady  wing  motion,  including  any  accelera¬ 
tion  of  the  wing  vertex  ae  well  as  pitch,  yaw  and  roll  about 
•pacified  axes.  Any  built-in  camber  of  the  wing  can  also  be 
included  in  the  form  (3.12).  The  angle  UK"  is  an  apparent 


“roll*  angle,  observed  at  fixed  z,  which  it  actually  a  com¬ 
posite  of  the  usual  Euler  angles.  For  this  reason  Xo,  Y0  and 
K  can  depend  on  *  as  wsll  as  t.  For  simplicity  we  assume 
the  wing  trace  to  be  fiat,  although  the  method  used  here  can 
be  generalised  to  include  any  moderately  thin  cross  section 
as  well  as  control  surface  deflections.  Note  that  j),  in  this 
case,  is  the  instantaneous  spanwise  coordinate  as  shown  in 
Figure  7,  and  £  is  the  effective  semispan  at  that  instant  at 
the  specified  z.  Thus,  Xo(z,t),  Yo(z,t)  give  the  instanta¬ 
neous  position  of  the  mid-effective  span  which  is  not  on  the 
wing  center  line  during  yaw. 

If  there  were  no  wake  the  discontinuous  part  of  the  ve¬ 
locity  at  the  wing  trace,  AI?a,  would  be 

A^oltrace  = 

20  r /„  DX0\_„  „D.Y0} 


«¥  m-r- 


r  Irr  .  —  DXo\  __  .  DMYo 1 
^C/oo  sin  a - cos  K  -  smK  — -I 

l>2  ,,  ($)  \ 


t  -p 


since  in  that  case,  according  to  (3.10),  we  would  have 
Trsxcs  =  0.  For  example,  in  the  limit  of  steady  motion  at 
small  a  and  with  no  yaw,  pitch  or  roll,  (3.13)  reduces  to 


^oltrace  —  - ' 


T~y2 


(3.13a) 


corresponding  to  Jones’  classical  result  [11],  Note  that  in 
either  (3.13)  or  (3.13a)  the  flow  predicted  in  this  way  (with¬ 
out  vortical  wakes  above  the  wing)  is  singular  at  the  wing 
edges.  The  next  step  is  to  add  the  wakes  and  to  determine 
their  vorticity  distributions  so  as  to  ‘‘smooth’’  the  flow  at 
the  edges. 

In  a  manner  similar  to  that  used  for  the  2D  airfoil,  we 
introduce  an  “angle  variable*  such  that 


y  = -sin/3  (3.14) 

and  use  the  theory  of  conjugate  functions  to  write  down 
the  additional  vorticity  (or  jump  in  1^)  required  at  the  wing 
trace  in  order  to  continue  to  satisfy  the  boundary  condition 
(3.12)  despite  the  velocity  field  (3.11)  induced  by  the  wakes. 
The  result  is 

A^ltrace  =  A^oltrace  +  ‘^iltrace 


‘urace  -  ''  ? - - — 

to  cos  0 

,  2tan/3 1*  ^Mwake  co«  r) 

T  COS  T  —  COS  0 


The  wing-edge  singularities  in  these  additional  “wake- 
induced*  terms  must  exactly  cancel  the  singularities  in 


APo  Itracs  at  both  edges  of  the  wing  if  our  Kutta-condition 
postulate  is  to  be  satisfied.  Note  that  at  this  point  we 
can  no  longer  simply  assume  rrexcs  =  0,  but  must  apply 
(3.10)  in  order  to  try  to  determine  it. 

In  Eqn.  (3.15)  Qwabe(P,  t)  can  be  written  down  in  terms 
of  the  wake  integrals  in  (3.11).  Cancellation  of  the  singu¬ 
larities  at  both  wing  edges  then  requires 

for  p  =  +f: 

,  [7..  .  _  DX 0\  „  .  „D,Fol  j D,K 

2  [(Resina  -  —  J  cos  K  -  sm K  —  j  +6-^ 
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fory  =  -f: 


_  2rrHAcm{*t  t) 
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tfFORT 
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f{2:°*r)*  +  %  y{trao)2  +  \ b 


In  the  above  relatione 


zrRr  =  {zrKT(u,t,z)  -  (Xo + iYo)) 


2r,o  s  ,-Jrfs.o (*««(,..  t|g)  -  (Xo  +  iYo)) 

with  the  Z.  as  defined  below  (3.11). 

Taking  the  difference  of  the  two  equations  involving 
Truces  yields  rT«>»c»  =  (I|~)  %r>  80  ‘hat  we  obtain 
the  constraint 

rro«r(*,0  +  rsrao(r,0  = - — - (3-1®) 

The  sum  of  the  two  relations  then  yields  an  equation  for 
the  vorticity  strengths  (circulation  elements)  in  the  wake 
traces  which  is  analogous  to  the  Wagner  integral  equation 
(3.9).  The  result  is 

(  tt  —  DXo\  —  tyD,Yo 
^1/ooSina - pj-)  C08^  -  sin  if  ^  — 


.  "  if*"*1 

=  ReZ.  /  iu  <  --  iu  ■  - +  — —  *-£■-■  } 

*Jo  +  J 

(3.17) 

The  left-hand-side  of  (3.17)  is  known  at  each  z  and  t 
for  any  given  imposed  wing  motion.  Thus,  (3.17),  together 
with  (3.16)  and  (2.20),  determines  the  circulation  elements 
throughout  the  wake  system  within  the  present  model.  So¬ 
lution  proceeds  in  a  manner  mathematically  quite  similar 
to  the  2D  airfoil  case. 
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4  Loads,  Forces  and  Moments 

4.1  2D  Airfoil  Case.  Large* Amplitude 

Unsteady  Response 


In  R«f.  [1  ]  the  force*  end  moment*  occurring  on  20  air¬ 
foil*  in  large-amplitude  unateady  motion  were  calculated 
using  the  pressures  over  the  airfoil  aurface  a*  predicted  by 
the  Bernoulli  equation  for  unateady  potential  flow.  In  reck¬ 
oning  the  lift,  a  ‘leading-edge  force*  wa*  included  in  the 
classical  manner,  except  that  no  linearization*  were  admit¬ 
ted.  Thi*  procedure  wa*  built  into  the  “NLWAKE*  code 
reported  in  [l|. 

More  recently,  we  have  reformulated  the  lift  and  moment 
expreaaion*  by  taking  advantage  of  the  possibility  of  ex¬ 
pressing  all  or  moat  wake-related  term*  (without  linearii- 
ing)  in  term*  of  the  corresponding  “wake  integral**  -  in¬ 
tegral*  over  the  (known)  wake  vorticity.  Thi*  itep  mimic* 
classical  procedure*  in  the  linear  theory  [18]  and  establishes 
the  equivalence  of  the  original  formulation  in  Ref.  [1]  with 
the  application  of  the  principle  of  conservation  of  impulse. 
Most  importantly,  the  reformulation  greatly  improve*  the 
numerical  accuracy  achievable  in  the  codes.  The  reason 
is  that  several  large  contributing  term*  in  Ct  and  Cw,  in 
XLWAKE  almost  exactly  cancel,  leaving  crucial  small  dif¬ 
ference*  which  can  be  difficult  to  calculate  with  sufficient 
accuracy.  In  the  extension  of  NLWAKE  to  include  the  new 
formulation,  these  partial  cancellations  are  explicit  and  the 
remainders  are  written  in  a  manner  more  suited  to  precise 
numerical  evaluation. 

The  (non-linear)  expressions  for  lift  and  moment  includ¬ 
ing  the  appropriate  wake  integrals  in  this  way  are  written 
out  here  both  for  completeness  and  for  comparison  with 
their  classical  linear  versions.  For  the  lift  coefficient,  using 
the  notation  of  Eqn.  (2.11),  (3.3)  and  (3.5),  we  find 


~  ^twake(*) 

where  o  =  2 i,/c  with  2,  defined  as  below  (3.9),  and  U  = 
i'oo  -r  htana.  The  small  non-linear  term  Ct  waj.e(0  is 


i  waksvO  - 

2  fi 

c(fl>  cos  a  j  ^  **  ^wake  (*<  0  (7o(*.0  +  7t  (*,()) . 


more  useful  her*  simply  to  not#  that  Ct  wajj#  vanishes  iden¬ 
tically  in  the  linear  limits  of  a  wake  extending  straight  be¬ 
hind  the  airfoil  since  it  involve*  &Wak*'  Also,  we  have  not 
yet  encountered  a  case,  even  in  the  non-linear  treatment,  for 
which  Ct  wijj#  ha*  been  important,  even  though  it  clearly 
can  be  significant  or  even  crucial,  for  a  maneuver  in  which 
the  airfoil  file*  “through,*  or  just  above  or  below,  its  own 
wake. 

The  fourth  term  on  the  rhs  of  (4.1)  also  vanishes  in  the 
linear  limit,  since  in  that  case  both  j  (57)  and  U  cos  a  ap¬ 
proach  C/00.  This  term  in  fact  illustrates  the  numerical  ad¬ 
vantage  of  the  wake-integral  formulation,  since  it  represents 
explicitly  the  partial  cancellation  between  two  otherwise 
quite  large  terms. 

The  first  three  terms  on  the  rh*  of  (4.1),  on  the  other 
hand,  approach  their  classical  linear  counterparts  exactly 
in  the  limit  of  small  amplitude  motion.  The  2nd  term  is 
often  called  the  ‘apparent  mass*  term,  and  the  3rd  the 
“wake  effect.*  Thus,  on  using  (4.1),  we  know  that  the  non¬ 
linear  theory  appropriately  reproduce*  the  linear  version, 
and  includes  it. 

The  corresponding  moment  coefficient,  C„0 ,  about  mid¬ 
chord  is  given  in  similar  fahsion  by 


„  U  cos  a  f  *  it  t  7o(i,f) 
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where 


C, 


wake 


(Os 


idi  “wake(i-t)(7o(i,f)  -  7i(i,0). 


(4.2a) 


Once  again,  the  final  two  terms  on  the  rhs  of  (4.2)  vanish  at 
small  amplitudes,  and  the  first  three  terms  reduce  precisely 
to  their  classical  counterparts  in  the  linear  limit. 

Computer-generated  results  corresponding  to  Eqns.  (4.1) 
and  (4.2)  are  illustrated  in  Section  5  in  terms  of  2D  airfoil 
response  to  large-amplitude  imposed  motions. 


4.2  3D  Case.  Loads  and  Normal-Force 
Distribution  on  Slender  Wings 


(Via) 

and  can  also  be  written  out  in  terms  of  wake  integral*.  The 
resulting  expression  is  cumbersome,  however.  It  is  perhaps 


The  pressure  difference,  or  loading,  across  the  wing  for 
the  3D  case  can  again  be  calculated  using  the  Bernoulli 
equation  for  unsteady  potential  Sow  on  either  side  of  the 


wing.  A*  before,  the  wake- related  term*  are  two-fold.  First, 
the  aeeociated  velocity  field  of  the  wake,  ae  given  by  Eqn. 

(3.11) ,  contribatee  directly  to  the  preeeure  variation*.  Sec¬ 
ondly,  the  apparent  upwa*h  caused  by  the  wake  must  be 
cancelled  out  by  the  effect*  of  an  addition  wing-trac*  vor- 
ticity,  AVi(g,t),  a*  given  in  (3. IS),  thus  implying  an  ad¬ 
ditional  jump  in  the  term  in  the  Bernoulli  relation.  In 
calculating  the  latter  we  take  advantage  of  the  possibility 
of  writing  *(z,y,t)  on  the  wing  on  either  side  (+,  -)  as 

<fi±(z,y,t,t)  =  *(z,t;$,4  =0*) 

where  the  coordinates  (i,  y)  are  as  illustrated  in  Figure  7. 
We  find,  for  slender  wings,  with  V  =  (d,  0,  w) 

1  1  lt0P 

-Ap=  -p  =p“(z,t;y,0")  -p+(z,rj,0+) 

Po°  *  bottom 

-<K>  iK+i(Jt+»|i)  (4.3) 

a< V>AV 

where  =  £  +  W  -fa  as  before.  In  writing  (4.3)  we  have 
also  used  the  notation  defined  in  (2.6)  and  (2.7)  and  below 

(3.12) ,  as  well  as  in  Figure  7. 

As  mentioned  in  Section  (3.2)  proper  use  of  the  slender- 
wing  concept  allows  us  to  repace  VR  with  W ,  provided  we 
recall  that  the  low  /R  wing  is  “slender”  along  the  chordwise 
direction  and  not  necessarily  along  the  z-direction  during 
violent  pitch  and  yaw. 

The  A#  occurring  in  (4.3)  can  be  written 


DtX o  DtYo  zs j 


+C,wake  +  Cn 


C-wak.('-0  S  iklft  {lJ_\  ^ak. 
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J*d0 cos/??wake  (&v0  +  A?,)  (4.5 


Using  (3.11)  we  can  express  C„Wake(*>0  ‘n  t«rms  of  a 
wake  integral  analogous  to  the  2D  airfoil  case.  The  result 


=  ikwtRel 


a.  =  2 Zra°lb  I  =  2Z:°*T/b 


A<&(y,z,t)  =  rjrao(z,f) 


-J’dy  A V(5,z,i 


where  AV'  is  given  in  (3. *3)  and  (3.15).  Further,  <  ^  >  = 
-Uoo  sin  5  sin  K  •+•  Vwake. 

Expression  (4.3),  with  (4.4),  lend*  itself  to  relatively 
straightforward  calculation  of  the  normal  force  distribution 
C, (z,t).  Spanwise  integration  with  respect  to  $  proceeds 
in  a  manner  mathematically  the  same  as  for  the  2D  airfoil 
case,  except  that  in  the  present  model  no  edge  forces  are 
included.  The  lift  and  drag  and  also  the  roll,  pitch,  and 
yaw  moment*  of  a  given  wing  can  be  reckoned  in  a  simi¬ 
lar  manner,  on  including  subsequent  weighted  integrations 
ever  z.  We  find 


C„(z,t)  = 
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and  we  note  a,  =  -i  at  the  starboard  edge  and  a,  =  -t-t  at 
the  port  edge. 

In  the  interactive  method  for  the  slender  wing  case  this 
latter  “wake-integral"  is  easy  to  calculate  on  the  computer, 
since  the  necessary  information  becomes  available  automat¬ 
ically  at  the  same  time  the  wake  structure  itself  is  being 
generated  using  the  method  described  in  Section  (3). 

The  C„  ( z,  t)  term  can  also  be  calculated  in  terms  of  wake 
integrals,  and  that  approach  appears  to  be  the  most  con¬ 
venient  and  numerically  accurate  method  to  determine  its 
value.  We  note  here  only  that  C„  has  no  linear  counterpart 
(and  vanishes  by  symmetry  in  certain  simple  flight  config¬ 
urations),  yet  its  effect  appears  to  be  potentially  of  great 
importance  in  certain  maneuvers  involving  intimate  inter¬ 
action  between  the  wing  and  its  wakes. 

The  first  term  on  the  rhs  of  (4.5)  corresponds  for  small  3 
to  the  linear  theory  result  of  Ref.  fllj.  In  fact,  for  pointed 
slender  wings  in  steady  flight  at  small  3,  this  term  inte¬ 
grates  chordwise  to  yield 

Ct  =  ~ci,  (4.7) 

the  classical  Jones  result.  The  additional  terms  in  the  corre¬ 
sponding  curly  brackets  represent  the  local,  instantaneous 


normal-fores  effect*  of  plunge,  pitch,  etc.,  that  would  ariss 
if  there  were  no  wait**. 

The  wake-intsgral  term*,  40(1  <?•»(*, *)•  r#P" 

resent  “history  effect**  in  the  relatiomhip  between  wing 
motion  and  wing  loading  or  reepon**.  A*  we  have  *een, 
the  variou*  vorticity  element*  in  each  wake  reflect  condi¬ 
tion*  imposed  at  upstream  wing  edge*,  each  at  an  earlier 
time.  Thus,  for  example,  in  the  event  of  a  *udden  maneu¬ 
ver  of  short  duration,  the  outer  parte  of  the  wake  affecting 
<?„(*,()  will  often  represent  wing  attitude  and  other  condi¬ 
tions  in  effect  before  the  maneuver  began.  Conversely,  af¬ 
ter  a  given  maneuver  is  ‘over*,  many  outer  wake  element* 
will  still  reflect  the  actions  taken  during  the  maneuver,  so 
that  wing  loading  will  take  time  to  adjust  and/or  return  to 
“normal.*  Naturally,  any  occurrence  of  vortex  breakup  can 
have  a  drastic  effect  on  this  history- related  part  of  the  wing 
response,  although  it  may  not  always  be  unfavorable.  Of 
course,  before  we  can  regard  wing  response  a*  “determinis¬ 
tic*  we  must  learn  how  to  predict  both  the  occurrence  and 
specific  nature  of  such  vortex  bursting.  (12|-(14| 

5  Typical  On-Line  Results;  Non- 

Linear  2D  Airfoil  Response 

A  critical  part  of  studying  the  behavior  of  wings  in  ma¬ 
neuver  is  the  understanding  of  the  wing’s  aerodynamic  re¬ 
sponse  to  new  conditions  suddenly  imposed.  Representa¬ 
tive  of  such  response  are,  for  the  2D  case,  the  results  of  the 
classical  Wagner  problem,  illustrated  for  the  low-amplitude 
linear)  caae  in  Figure  8(a).  The  net  bound  circulation  on 
;he  airfoil,  To  •+■  Tj,  adjusts  only  gradually  to  the  suddenly 
imposed  change  in  angle-of-attack  as  repreeented  by  the 
jump  in  To-  The  lift  coefficient  jumps  to  i  its  eventual 
value  and  gradually  adjust*.  Note  that  the  relaxation  to 
•.he  quasi-steady  result,  is  very  slow.  On  the  Figures,  time 
:$  normalised  to  (  j)  {c/Ua <,)• 

The  corresponding  nonlinear  (large-amplitude)  case  is  il- 
1-istrated  in  Figure  9(a).  Note  that  the  response  is  very 

i.m.lar  to  the  linearized  result,  despit*  ths  relatively  severe 
i?f:rmation  of  the  wake  as  illustrated  in  Figure  9(b).  This 
appears  to  be  a  result  of  the  fact  that  the  wing  in  each 
:ase  is  adjusting  to  shed  vorticity  which  is  rather  strongly 
::r.centrated  toward  the  far  reaches  of  the  wake,  ultimately 
several  chord  lengths  downstream. 

Illustration  of  the  airfoil's  aerodynamic  response  to  a 
rudimentary  but  sever*  maneuver  is  illustrated  in  Figure 
There,  a  sudden  single-cycle  sinusoidal  angle  of  attack 
variation  of  large  amplitude  is  imposed  and  just  as  suddenly 
dipped.  The  wake  vorticity  coagulate*  and  begins  roll-up 
without  significant  net  displacement  from  the  plane  of  th* 
-...-foil  (Figure  10(a)).  The  lift  jumpe  quickly,  and  then  lead* 
the  quasi- steady  value  into  ths  negative  lift  region  (Figure 
r.'b)).  It  then  overshoote  the  recovery,  plunges  again,  and 
r..y  very  slowly  recovers  to  the  state  before  the  maneuvers, 
.i-strating  a  “history  effect*  for  th*  2D  cast  suggestive  of 
•he  phenomena  discussed  at  the  end  of  Section  4.  During 


the  maneuver,  from  start  to  finish,  the  airfoil  advance*  S 
chord  lengths  in  this  example.  Both  “apparent  mate*  ef¬ 
fects  and  upwash  associated  with  the  developing  wake  play 
vital  roles  in  the  net  reeult. 

The  pitching  moment  for  this  example,  though  large  in 
amplitude,  ii  eeientiaily  classical  in  response,  except  for  the 
relatively  pronounced  hietory  effect  a'ter  the  maneuver  is 
complete  (Figure  10(c)). 

A  number  of  additional  example#  of  airfoil  response  to 
imposed  large-amplitude  motion*  of  various  types  are  re¬ 
ported  in  Ref.  [l).  It  is  especially  instructive,  in  addition 
to  the  figure#,  to  actually  run  the  codes  of  (l)  on-line,  tak¬ 
ing  advantage  of  the  almost  instantaneous  wake  displays  to 
develop  a  sense  of  airfoil  response  and  the  reason*  for  it. 
Several  of  these  “flight*  sessions  have  been  taped  recently. 
In  the  3D  case  one  can  carry  out  similar  on-line  runs  with 
split-screen  displays  of  th*  wing  configuration,  wake  defor¬ 
mations  at  various  locations,  and  wing  loading  response. 
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1.  Particle  Path*  v*.  Streakline.  Formation  of  a  Vortex  Sheet. 

2.  Wake  Behind  an  OeciUating  Airfoil. 

3.  Sketch  of  Wakes  from  a  Delta  Wing. 

4.  Typical  Wake  Traces,  a)  Non-Linear,  b)  Linear,  Idealized 

5.  Wing  in  Maneuver,  a)  Wing  and  Wake  Trace*  at  Various  Chordwite  Locations,  b)  A  Streak  Surface,  Starboard 
Side,  for  given  t. 

6.  Superposition  Algorithm  for  Determining  Bound  Vorticity  Necessary  to  Cancel  Upwash  at  the  Plate. 

7.  Choice  of  Cross-Flow  Observation  Plane*,  a)  z  ~  z  Projection,  b)  x  -  y  Projections  and  “Roll*. 

8.  Airfoil  Response  to  Low-amplitude  Step  Change  in  Attitude,  a)  Circulation  and  Ct  vs.  Time,  b)  Wake  Display. 

9.  Airfoil  Response  to  Large  Amplitude  Step  Change  in  Attitude,  a)  Circulation  and  Ct  vs.  Time,  b)  Wake  Display. 

10.  Airfoil  Response  to  a  Single  Full-cycle  Completed  Imposed  Oscillation,  a)  Wake  Display,  b)  Lift  and  Circulation 
vs.  Quasi-steady  Values,  c)  Moment  Coefficient,  with  History  Effect. 
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Abstract 

An  aerodynamic  theory  of  the  flow  about  ilender  delta 
wing*  if  described.  The  theory  include*  a  treatment  of 
the  eelf-coneietent  development  of  the  vortex  wake  pattern* 
above  the  wing  neceeeary  to  maintain  emooth  flow  at  the 
wing  edge*.  The  paper  focuee*  especially  on  the  formation 
within  the  wake  of  vortex  ‘cores1*  as  embedded  regions  of 
total  pressure  lose,  fed  and  maintained  by  umbilical  vortex 
sheet*  emanating  from  the  wing  edgee.  Criteria  are  devel' 
oped  for  determining  the  growing  sise  and  location  of  these 
corse,  as  well  a*  the  distribution  and  strength  of  the  vortic> 
ity  within  them.  In  this  paper,  however,  the  poeeibility  of 
vortex  breakup  is  omitted.  The  aerodynamic  consequence* 
of  the  presence  and  evolution  of  the  core*  and  the  asso¬ 
ciated  wake  structure  are  illustrated  and  discussed.  It  is 
noted  that  wake  history  effects  can  have  substantial  influ¬ 
ence  on  the  distribution  of  normal  force  on  the  wing  as  well 
as  on  its  magnitude. 

1  Introduction 

The  underlying  assumption  adopted  in  this  work  is  that 
at  the  Reynolds  numbers  of  full-ecale  flight  it  should  be 
possible  to  advance  our  understanding  of  the  aerodynamic 
behavior  of  wings,  even  in  severe  maneuver,  without  nec¬ 
essarily  resorting  to  the  full  Navier-Stokes  equations.  This 
classical  notion  envisions  the  flow  past  the  wing  as  con¬ 
sisting  of  large  surrounding  ‘outer*  region*  of  essentially 
inviscid  flow  in  which  are  embedded  "inner*  flows,  many 
key  features  of  which  are  viscous  or  viscous-related.  These 
embedded  flows  may  include  not  only  the  usual  boundary 
layer*  and  viscous  wakes  but  also  larger  recirculating  re¬ 
gions  (such  as  separation  bubbles)  as  well  as  shed  vortices 
and  vortex  ‘streets*  complete  with  their  own  inner  ‘cores*. 
The  location,  sise,  and  shape  of  the  inner  embedded  flow 
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regions  is  determined  interactively  with  the  larger-scale  in¬ 
viscid  flow,  but  their  very  existence  depends  on  the  pres¬ 
ence  of  viscous  action  occurring  somewhere,  especially  at 
the  various  interfaces  with  the  outer  flow. 

To  the  extent  that  such  an  hypothesis  is  viable  one 
can  contemplate  an  analytic  approach,  albeit  computer- 
assisted,  to  the  study  of  non-linear,  unsteady  wing  aerody¬ 
namics.  In  many  situations  the  outer  flow  will  be  potential; 
even  the  embedded  flow  regions,  when  of  large  enough  scale, 
may  be  eeeentially  inviscid  within,  although  they  must  then 
be  surrounded  by  viscous  layers  or  Sheaths*  separating 
them  from  the  outer  (potential)  flow.  A  common  feature 
shared  among  the  embedded  flow  regions,  however,  is  that 
they  will  not  be  irrotational;  reduced  total  pressure,  dis¬ 
tributed  variously  over  the  inner  regions,  will  be  present. 
In  addition,  when  any  of  the  embedded  flow*  detach  from 
the  wing  surface  a  key  part  of  the  aerodynamic  problem  is 
to  determine  their  evolution,  in  sise,  shape  and  location,  as 
affected  both  by  convection  in  and  viscous  interaction  with 
the  surrounding  flows.  When  compressibility  is  important, 
the  possible  formation  of  shocks  as  well  as  the  propagation 
of  acoustic  disturbances  can  also  affect  the  interplay  of  the 
inner  and  outer  Sows.  In  this  paper,  however,  we  limit 
ourselves  to  th*  incompressible  case. 

If  flow  models  of  wing  aerodynamics  are  constructed 
along  these  lines,  fully  three-dimensional  analytic  studies  of 
the  problem  can  be  carried  out  and  the  procedure  is  essen¬ 
tially  classical.  Given  a  wing’s  geometry  and  attitude  and 
its  flight  statu*  or  sequence,  the  precise  boundary  condi¬ 
tions  needed  to  determine  the  flow  can  be  written  down  and 
classical  potential  flow  techniques  applied.  The  inner  flows 
can  be  treated  using  a  variety  of  approximate  treatments. 
The  algebraic  complexities  in  the  3D  case  are  formidable, 
however,  and  (at  least  at  the  outset)  may  serve  only  to 
obecure  the  progrees  we  wish  to  make  with  an  analytic  or 
quasi-analytic  approach.  In  this  paper,  therefore,  we  adopt 
the  idea*  of  the  slender  wing  theory  of  R.T.  Jones  [l]  as 
a  mean*  of  illustrating  the  proposed  techniques  in  a  some¬ 
what  simplified  framework.  Jones’  qua*i-2D  treatment  for 
low  aspect-ratio  wings  is  a  remarkably  powerful  and  insight¬ 
ful  tool  which  historically  ha*  also  provided  a  framework  for 
understanding  the  results  and  trend*  of  more  exact  numer¬ 
ical  studies. 

Advance*  beyond  Jones’  theory,  especially  a*  applied  to 
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delta  wing*,  began  with  the  recognition  by  many  authors 
[2],  [3|,  [4],  [S]  of  the  ne«d  for  a  vortex  wake  structure 
above  the  suction  surface  of  the  wing.  (Fig.  1)  The  role 
of  this  wake  pattern,  somewhat  idealised  here  in  the  fig¬ 
ure,  is  to  provide  means  of  establishing  smoother  flow  (a 
“Kutta  condition*)  at  the  wing  edges  than  that  arising  in 
the  original  Jones  treatment.  The  resulting  modifications 
provide  predictions  more  commensurate  with  real  flows  at 
large  but  finite  Reynolds  Nr.  One  can  view  such  a  dou¬ 
ble  wake  structure  as  the  result  of  a  form  of  leading-edge 
separation  occurring  along  each  of  the  swept  edges. 

In  this  paper  we  discuss  the  self-consistent  development 
of  the  wake  pattern  above  a  low  aspect-ratio  delta  wing  at 
modest  angle  of  attack  in  the  context  of  its  effect  on  slender 
wing  theory.  In  particular,  we  undertake  to  include  in  our 
treatment  the  presence  of  embedded  cores  of  distributed 
vorticity,  or  regions  of  reduced  total  pressure,  near  the  end 
of  each  wake  trace  as  illustrated  in  Fig.  2.  The  physical 
reasons  for  the  presence  of  these  embedded  regions  are  dis- 
cuseed  briefly  in  the  next  Section.  We  will  find  in  later  Sec¬ 
tions  that  the  effect  of  embedded  low-total-pressure  regions 
on  the  development  along  the  wing  of  various  flow  quanti¬ 
ties,  including  static  pressure,  requires  careful  attention.  A 
self-consistent  procedure  for  determining  these  parameters 
in  the  slender-wing  Limit  is  developed  in  this  paper  and 
compared  with  the  classical  case  of  purely  potential  flow. 
In  addition,  we  derive  criteria  for  determining  the  site  and 
evolution  of  the  vortex  cores.  However,  in  this  paper  we  do 
not  consider  the  possibility  of  vortex  breakup.  Some  of  our 
recent  work  on  that  subject  will  be  reported  elsewhere. 


2  Wake  Formation  in  Potential  Flow 

2.1  Development  of  Wakes  with  Total 
Pressure  Loss 

In  Reference  [6]  the  self-consistent  development  of  vor¬ 
tex  patterns  above  slender  delta  wings  in  large- amplitude 
unsteady  motion  is  discussed  from  the  point  of  view  of  an¬ 
alysing  the  cross-flow  pattern  at  each  chordwise  station.  In 
that  work,  however,  the  wakes  are  treated  ideally,  as  if  they 
were  simply  vortex  sheets  of  complex  shape,  with  no  total 
pressure  loss  in  the  flow.  The  focus  there  is  on  determining 
interactively  the  shapes  of  the  two  required  wakes  as  viewed 
in  successive  cross-flow  planes.  The  distribution  of  vortic¬ 
ity  along  them  is  determined  so  as  to  maintain  smooth  flow 
along  the  wing  leading  edges.  The  method  extends  the  clas¬ 
sical  literature  (2),  [3],  [4],  (5|  so  as  to  allow  for  severe  wing 
motion  during  maneuver.  Non-linear  self-consistent  vortex 
wake  evolution  is  explicitly  included. 

As  pointed  out  in  that  paper,  such  a  model  of  the  wake 
structure  is  somewhat  too  idealised,  in  that  the  develop¬ 
ment  during  wake  roll  up  of  vortex  cores  at  finite  Reynolds 
numbers,  with  attendant  regions  of  distributed  total  pres¬ 
sure  loss,  is  not  addressed.  Similarly,  the  related  occur¬ 
rences  of  secondary  separation  'bubbles,"  which  arise  due 


to  localised  separation  of  the  wing  boundary  layers,  and 
the  phenomena  of  vortex  breakup,  are  explicitly  omitted. 

In  the  course  of  the  non-linear  evolution  of  the  (ideal, 
inviscid)  vortex-sheet  patterns  as  in  [6j,  the  wakes  gener¬ 
ally  tend  to  "coagulate,"  with  the  vortex  traces  spiralling 
in  on  themselves  and  "rolling  up."  (See  Fig.  2)  If  this  is 
the  case  and  the  Reynolds  number  is  finite  (no  matter  how 
large),  viscous  dissipation  and  smearing  out  of  the  vorticity 
must  become  important  at  some  stage  as  the  contortions 
of  the  ideal  wake  traces  become  more  and  more  intricate. 
The  result,  as  sketched  in  Fig.  2,  is  the  creation  of  a  re¬ 
gion  of  smeared-out  vorticity  in  the  "core*  of  each  wake. 
(Multiple  cores  can  also  occur.)  In  these  core  regions,  a 
corresponding  embedded  region  of  total  pressure  deficit  is 
created.  Both  numerical  and  experimental  evidence  exists 
for  such  embedded  regions;  see,  for  example,  references  [7], 

[8!.  [91.  [10],  (llj. 
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Although  the  wakes  may  partially  roll  up  into  "cores* 
as  discussed  above,  we  still  require  a  careful  treatment  of 
those  portions  of  the  wake  which  remain  thin  and  sheet-Uke, 
with  negligible  total  pressure  loss  across  them.  For  exam¬ 
ple,  the  strength  of  the  cores  must  change  as  ‘.hey  move 
downstream  over  the  wing,  not  least  because  of  chordwise 
variations  in  wing  geometry.  Thus,  sheet-like  portions  of 
each  wake  can  be  regarded  as  "umbilicals*  feeding  each 
core  (see  Figures),  and  to  these  the  analysis  of  [6j  applies 
almost  without  change. 

In  the  present  paper,  however,  we  wish  to  take  advan¬ 
tage  of  a  form  of  the  time-analogy  approach  inherent  in  the 
classical  Jones-Munk  treatment  of  slender  wings  and  bod¬ 
ies.  [  1 ) ,  [  12]  For  simplicity  in  the  present  discussion  we  take 
the  wing  to  be  a  delta  of  planar  geometry  and  in  steady  mo¬ 
tion  at  small,  constant  angle  of  attack,  omitting  the  radical 
maneuvers  discussed  in  [6]. 

Jones  suggested  an  observation  coordinate  frame  fixed  in 
the  still  atmosphere  through  which  the  wing  being  observed 
penetrates  during  its  flight.  Extending  the  ideas  of  Munk, 
Jones  pointed  out  that  if  the  wing  was  slender  enough  (i.e., 
of  low  aspect  ratio)  an  observer  fixed  in  this  frame  would 
see,  in  a  “cross-flow  plane*  oriented  normal  to  wing  surface, 
a  quasi-2D,  but  unsteady,  flow  pattern.  The  unsteadiness 
would  arise,  for  example,  because  the  observer  would  see, 
in  the  cross-flow  plane,  a  wing  trace  "growing*  in  span  ac¬ 
cording  to 

6(f)  =  =  Uoo  COS  w  Uao  (2.1) 

■  2  2 

In  addition,  in  Jones’  scheme,  the  wing  would  appear  to  be 
plunging  at  constant  rate  Uoo  sin  a  ks  U^a. 

As  illustrated  in  Fig.  3,  we  adopt  in  this  paper  the  Jones 
"time-analogy*  view,  with  the  exception  that  we  "stop* 
the  apparent  plunging  by  including  a  countering  croes-flow 
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We  make  uee  of  Eq.  (2.3)  in  writing  out  the  “wake- 
history*  integral*  below.  In  [6j  it  ie  ihown  that  A  it  a 
“purely-convected*  quantity,  acting  at  a  label  for  a  given 
circulation  element  in  the  wake. 

2.4  Determining  the  Vortex  Strength  Dis¬ 
tribution  in  Ideal  Wakes 

The  net  circulation  in  the  crow  flow  plane  is  zero  for  all 
timee  by  Kelvin’s  Theorem; 

psrtD  +  peoer  +  preset  _  0  (2 .4) 

where  rTKACa  is  the  crost-flow  circulation  of  the  wing 
trace,  and  r*r*°,  r,°',T  refer  to  the  starboard  and  port 
wake  traces.  From  the  symmetry  of  the  problem  with  zero 
roll  angle  we  see  that 


When  the  wing  it  in  steady  flight  at  angle  of  attack  a,  we 
can  relate  time  variations  seen  in  the  present  observation 
frame  to  chordwise  variations  referred  to  in  [6]  through  the 
transform 

h=U~e° taTzt*U-J-z  (2-2) 

2.3  Convection  of  Ideal  Wakes 

The  mathematical  description  of  the  non-linear  convec¬ 
tion  and  evolution  of  ideal  wakes  it  contained  in  Reference 
|6|  and  will  not  be  repeated  here.  However,  a  synopsis  of 
the  general  approach  it  uteful. 

In  formulating  the  ideal  wake  model  the  only  concession 
to  the  existence  of  viscous-related  phenomena  is  to  enforce 
a  Kutta  condition  •  i.e.,  smooth  flow  -  at  each  edge  of  the 
wing  trace.  Once  an  element  of  circulation  hat  been  shed,  in 
that  model,  the  wake  it  treated  in  the  inviscid  limit.  Thus 
each  increment  of  vorticity  it  convected  in  accordance  with 
the  inviscid  Helmholtz  relation. 

Moreover,  the  ideal  treatment  requires  that  each  wake 
trace  be  ■force-free’  along  its  entire  length.  The  wake  vor¬ 
ticity  then  convects  at  the  local  flow  velocity.  In  the  cross¬ 
flow  plane  this  convection  velocity  at  any  point  on  the  wakes 
is  the  sum  of  the  crose-flow  component  of  the  free  stream 
velocity,  t/«>  sin  a,  and  the  "induced*  velocities  associated 
with  all  the  neighboring  vorticity  on  the  wing  trace  and  in 
the  port  and  starboard  wakes. 

Owing  to  the  wake  history,  each  element  of  a  specified 
wake  trace  (i.e.,  port  or  starboard)  can  be  identified  accord¬ 
ing  to  the  time  it  emanated  from  the  wing  edge.  Denoting 
this  time  at  A,  the  two  wake  trace*  at  each  instant  t  are 
streaklines  along  which  A  varies  from  0  at  the  far  reaches 
of  the  wake  trace  to  t  at  the  wing  trace  ends.  (See  Fig.  4) 
Hence  the  incremental  circulation,  dT,  can  be  written 

dT=^dX.  (2.3) 


^wake  ’^wake  —  ^wake(^) 

'  iH  {/'"  (!  6) 

In  writing  Eq.  (2.6)  we  have  made  use  of  Eq.  (2.3).  In 
addition,  ZBOKT{ A)  and  Z?TBD( A)  are  the  instantaneous 
complex  locations  in  the  cross  flow  plane  (x,  y)  of  the  vor¬ 
ticity  elements  in  the  port  and  starboard  wake  traces,  re¬ 
spectively. 

For  our  chosen  reference  frame,  the  boundary  condition 
of  no  flow  through  the  wing  trace  is  written 

^Iwing  trace  =0>  “  J  -  *  -  J  '  1  =  0  '  (2-7) 

In  Eq.  (2.7)  y  is  the  instantaneous  spanwiee  coordinate,  and 
|  is  the  semispan  at  the  specified  t.  (Fig.  4). 

We  define  the  complex  velocity  Wa  =  Uo  -  »V0  to  be  that 
which  would  occur  if  there  were  no  wakes.  In  that  cate, 
imposition  of  boundary  condition  Eq.  (2.7)  in  the  presence 
of  the  cross  flow  I/*  sin  a  leads  to  a  discontinuous  part  of 
the  velocity  at  the  wing  trace,  A  Vo: 


=  -2Uoo  tin  a  - 


This  it  Jones’  classical  result  in  the  limit  of  small  angle  of 
attack,  sin  a  or  a.  It  is  alto  singular  at  the  wing  edges. 


psrto _  _rfo*f 

p TKACM  _  Q  (2.5) 

(See  Fig.  4)  (For  cates  involving  asymmetric  maneuver 
preset  mugt  worked  out,  and  will  not  generally  vanish.) 

In  complex-variable  notation  the  quati-2D  cross-flow  ve¬ 
locity  associated  with  the  wakes  it,  with  Z  =  x  +  iy 


velocity  (from  the  left,  in  Fig.  3)  of  magnitude  ffootina. 
Thus,  in  th*  cross-flow  observation  plane,  the  origin  of  the 
x  -  y  frame  remain*  fixed  on  th*  wing  center  line  at  the 
wing  passes  throagh. 

It  it  in  describing  th*  evolution  of  the  two  wake  traces 
that  the  time  analogy  approach  ha*  its  greatest  appeal.  In 
the  spirit  of  Jones’  low  aspect-ratio-limit  treatment,  the  po¬ 
tential  flow  solution  in  the  observation  plane  at  each  instant 
represents  an  instantaneous  solution  of  Laplace’s  Equation 
corresponding  to  the  boundary  conditions  of  the  moment. 
But  these  now  include  not  only  the  apparently  growing  wing 
trace,  but  also  the  developing  wake*  (or  wake  traces)  cate¬ 
nating  from  the  wing  edge*.  Moreover,  at  any  instant  t  the 
outer  portion*  of  each  wake  trace,  including  and  approach¬ 
ing  the  cores,  carry  with  them  the  “history*  of  the  wing’s 
penetration  of  our  observation  plane;  the  vorticity  in  those 
parts  of  the  wake  were  generated  at  times  previous  to  t. 
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Next,  we  introduce  the  complex  velocity  3  Ui  -  i  Vt 
to  represent  the  how  held  associated  with  the  additional 
vorticity  at  the  wing  trace  needed  to  continue  to  eatiefy 
Eq.  (2.7)  in  the  pretence  of  ae  given  by  Eq.  (2.6). 

Ueing  the  theory  of  conjugate  functione  we  And 


^^‘Itrace 


2  tan  0  1*  dr(-Uwike  coe  r) 
x  J  _  ^  coe  r  -  cos  0 


(2.9) 


where  AK,  it  the  jump  in  Vi  at  the  wing  trace.  In  Eq. 
(2.10)  the  “angle  variable*  0  it 


y=b-*\n0  (2.10) 

with  0  being  §  and  -f  at  the  port  and  starboard  wing 
edges,  respectively. 

The  wing-edge  singularities  in  these  additional  ‘wake- 
induced*  terms  must  exactly  cancel  the  singularities  as 
given  in  Eq.  (2.8)  at  both  edges  of  the  wing  if  the  Kutta- 
conditions  are  to  be  satisfied. 

In  Eq.  (2.9)  UW3^t  can  be  written  down  in  terms  of  the 
wake  integrals  in  (2.6).  Cancellation  of  the  singularities  at 
both  wing  edges  (y  =  +j,  -  j)  then  requires 


This  equation  provides  us  with  the  means  to  calculate 
the  rate  at  which  vorticity  is  shed  from  the  wing  edges 
while  satisfying  both  Kutta  conditions.  Its  derivation  is 
analogous  to  that  used  by  Wagner  [13]  in  developing  an 
integral  equation  for  the  wake  vorticity  in  unsteady  airfoil 
problems.  For  this  reason  we  refer  to  Eq.  (2.11)  as  the 
“Cross  Flow  Wagner  Integral  Equation*. 

The  left  hand  side  of  Eq.  (2.11)  is  known,  being  sim¬ 
ply  the  croee-fiow  component  of  the  free  stream.1  Opera¬ 
tionally,  on  the  computer  the  use  of  Eq.  (2.11)  is  relatively 
straightforward,  except  for  finding  the  initial  incremental 
vortex  elements  shed  at  the  instant  the  point  of  the  delta 
pierces  the  cross  flow  plane  (i.e.,  at  t  =  0).  Difficulties 
arise  there  since  the  kernel  of  the  rh*  is  singular.  An  effec¬ 
tive  method  for  treating  this  critical  first  moment  can  b>‘ 
developed  in  analogy  with  the  method  given  in  [14].  After 
the  procedure  has  been  properly  started  it  proceeds  without 
further  difficulty.  Once  the  strength  of  a  vortex  element  has 
been  determined  by  Eq.  (2.11)  at  the  instant  it  leaves  the 

lTh«  IKt  of  Eq.  (2.11)  it  really  »  »tep  function,  being  tero 
for  t  <  0  and  C/<»eina  thereafter.  For  the  prerent  care,  if  a  ir 
limited  to  email  angler  of  attack.  Um  tin  a  »  t/» a. 


wing  trace,  that  strength,  as  measured  by  its  incremental 
circulation,  remains  unchanged  as  the  element  convects  in 
the  wake.  Throughout  the  evolution  of  the  wake  we  “track* 
the  locations  of  the  elements  by  the  method  of  [6]  and  de¬ 
velop  a  “history*  of  the  wake  vorticity  density.  Thus  all  the 
information  needed  to  carry  out  the  wake  integrations  (such 
as  in  Eq.  (2.11))  is  present  except  for  that  part  which  de¬ 
pends  on  the  elements  currently  being  shed  into  each  wake. 
Eq.  (2.11)  then  determines  these  additional  entries  so  that 
the  edge  Sow  remains  smooth. 

2.5  Extension  of  the  Cross-Flow  Wagner 
Integral  Equation  to  Flows  with  Dis¬ 
tributed  Vorticity 

Eq.  (2.11)  can  be  reformulated  to  incorporate  wakes  in 
which  there  are  regions  of  smeared  out  vorticity  as  well  as 
“ideal*  regions.  As  an  intermediate  step  in  this  reformula¬ 
tion  we  simply  write  the  expression  for  the  wake  “induced* 
velocity  components  in  terms  of  integrals  taken  over  the  in¬ 
stantaneous  lengths  of  the  umbilicals  (i.e.,  those  portions  of 
the  wake  which  are  thin  and  sheetlike)  and  over  the  areas 
of  the  vortical  cores.  (See  Fig.  5)  Along  the  umbilicals  the 
incremental  circulation  is  written 

dT=^dt  (2.12) 

where  it  is  the  differential  length  measured  along  the  in¬ 
stantaneous  trace.  In  the  cores  the  incremental  circulation 
is  given  by 

dr  =  n<M  (2.13) 

where  O  is  the  vorticity  density  and  4A  the  differential  area. 
(See  Fig.  5)  Recalling  that  the  Biot-Savart  law  applies  as 
well  to  distributed  vortical  regions  we  obtain  the  “induc¬ 
tion*  equation 


^wake  ~  ^wake  ’"'’wake  — 


where  “£umb"  and  “Acore*  represent  integrations  along 
the  umbilicals  and  over  the  areas  of  the  cores  respectively. 

Since  the  shapes  and  strengths  of  the  umbilicals  are  not 
known  o  priori  we  again  use  the  concept  of  history  integrals 
and  ideal  wake  tracking  and  obtain: 


^wake  —  ^wake  ’"''wake  — 


» 


The  lower  limit  of  the  history  integrals  A:unct  denotea 
the  circulation  element  of  the  enda  of  the  umbilicals  at  the 
juncture  with  the  coree  while  the  upper  limit  t  pertaina  to 
the  vorticity  currently  entering  the  wake  at  the  wing  edgee. 
The  value  of  Xjnnct  iacreaeee  in  time  a a  vorticity  from  the 
umbilicala  *feeid*  the  coree,  and  muat  be  determined.  (See 
§  5). 

Applying  the  boundary  condition  of  no  flow  through  the 
wing  trace  Eq.  (2.7),  and  enforcing  the  Kutta  conditions 
we  obtain  in  a  manner  analogoua  to  that  of  $  2.4: 


Uoo  sin  a  = 


I  rt  i 

_  [  *raTBD  ) 

*junct  1  J{Zr‘D?+ki\ 


f aFORT 

ACOT9 


.(  ;■/,  °'m 
J*cor\D  {  J{ZST “>)'  +  $ 


Denoting  the  net  circulation  in  the  core  aa  Tcore.  and  that 
in  the  umbilical  aa  rumb,  the  total  circulation  in  the  wake 
ia  written 

^wake  =  Tcore  +  Tumb  ■  (2.17) 

Aa  pointed  out  in  the  Introduction  the  formation  of  em¬ 
bedded  vortical  regions  ia  dependent  on  the  action  of  vis¬ 
cosity  no  matter  how  large  the  Reynolds  number.  In  those 
regions  the  ideal  wake  tracking  method  ia  no  longer  directly 
applicable.  Therefore  we  need  to  aise  and  locate  each  core 
and  determine  its  vorticity  distribution  before  the  extended 
integral  equation  Eq.  (2.16)  can  be  helpful  to  ua.  We  show 
in  Section  6  that  one  muat  consider  the  role  of  viscous  in¬ 
teraction  between  the  corea  and  surrounding  flow  and  also 
include  criteria  to  determine  the  rate  at  which  vorticity  ia 
removed  from  the  umbilical.  A  sample  calculation,  includ¬ 
ing  the  sising  of  the  corea,  ia  described  there. 


3  Slender  Wing  Expansion 

3.1  Ordering  in  the  Slenderness  Parame¬ 
ter:  Formal  Procedure. 

Suppose  we  flx  ourselves  in  a  reference  frame  defined  by 
£,,  £»>  £*I  where  £,  and  £,  are  the  unit  vectors  in  the 
cross-flow  plane  depicted  in  Fig.  3,  and  £«  —  £«  *  £»  >* 
aligned  with  the  center-line  of  the  wing.  In  this  reference 
frame  the  wing  appears  to  travel  in  the  -z-direction  at  a 
speed  cos  a.  In  addition  there  is  the  steady  cross-flow, 
Ua o  sin  a  described  above. 


We  write  the  velocity  vector  v,  in  the  form 
JC  =  SL-l  -t- 


JU  =  («■»)  =  «£.  +  «£»  (3-2) 

are  the  cartesian  velocity  components  in  the  cross-flow 
plane.  Similarly,  if  D  =  curl  v,  is  the  vorticity, 

a  =  n.£x  +  ax  (3.3) 


Vxxvx-n,£,_  (£ l~^)i ••  V^-i‘-h+i’dy 

ax-(n.,o.);n.-£-£;n,«£-|j  (3.4) 

We  next  introduce  the  slenderness  parameter,  e,  charac¬ 
terising  the  wing  aspect  ratio,  /R.  We  can  then  write 

|^|/|Vx|  =0(0.  (3.5) 

In  the  light  of  Eq.  (3.5),  the  order  of  magnitude  cl  vari¬ 
ous  terms  in  Eqa.  (3.1)  through  (3.4)  can  be  worked  out. 
Baaed  on  the  time-analogy  approach  discussed  in  Section 
2.2,  chordwise  variations  are  related  to  time  variations  as 
in  Eq.  (2.2).  Eq.  (3.5)  is  then  equivalent  to 

|^[/(C/oo|Vx|)  =  cos  a  0(e) 

=  0(e)  (3.6) 

since  Uoo  is  of  order  unity.  Note  also  that  div  O  =  0,  so 


and  consequently 


|nx|/|n,|  =0(e). 


Eqs.  (3.4),  (3.5)  and  (3.8)  then  imply 


S=ow- 


The  above  corresponds  to  the  standard  procedure  show¬ 
ing  how  slenderness  is  employed  to  make  approximations 
for  low  /R  wings;  the  technique  is  used  throughout  this 


3.2  Flow  Quantities  in  the  Cross-Flow 
Plane  Through  First  Order. 

§  3.1  establishes  a  formal  procedure  for  calculating  the 
magnitude  of  various  flow  quantities  under  the  slenderness 
approximation.  One  begins  by  considering  the  cross-flow 
plane.  In  this  section  we  make  use  of  slenderness  to  derive 
the  basic  governing  equations  in  that  plane,  including  first- 
order  "time’-variations  in  the  time-analogy. 
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The  continuity  equation  in  incompreeaible  flow  is 


% 


V*=0.  (3.10) 

Since  =  0(<3)  Eq.  (3.10)  reduce*  to 


du  9v 
9x  +  3jr 


S  Vx  £x  =0  . 


(3.11) 


through  first  order  in  c.  This  can  be  satisfied  at  any  instant 
using  a  streamfunction  rp(z ,  y;  f)  such  that 


»i=(t.,»,)  =  Viii>  x  «,  (3.12) 

to  the  required  order.  (Note  that  £x  =  Uaa  sin  a  £,  ss 
UaoCit',  far  away  from  the  wing  trace  in  the  cross- flow 
plane.)  The  e-component  of  the  vorticity  is  then,  to  the 
same  order, 

n„  = -Vi i(>  .  (3.13) 


The  unexpanded  momentum  equation  in  the  “inviscid" 
regions  described  in  the  Introduction  can  be  written  as 

^  +  «  V£  =  ~iVp.  (3.14) 

Again  using  slenderness,  Eq.  (3.14)  becomes,  through  first 
order, 

^  +  £x  •  Vx  £x  =  Vx  P  •  (3.15) 

Here,  "pm  is  the  (constant)  fluid  density  and  “p*  the  static 
pressure.  Eq.  (3.15)  can  also  be  written  as  [I5j 


(3.16) 


The  usefulness  of  this  equation  will  be  seen  in  $  5. 

Utilising  Eq.  (3.11),  we  also  have  V3  ^  =  0,  defining  the 
usual  potential  problem  in  the  cross-flow  plane,  [lj,  [12]. 

By  contrast,  in  the  “inner*  vortical  embedded  regions 
(see  \  1),  even  if  they  are  essentially  inviscid,  we  must  adopt 
a  somewhat  different  approach.  Using  Eq.  (3.12),  we  can 
rewrite  Eq.  (3.17)  in  this  case  as 

-  iVxPx „  =  n.Vxtf  +  Vx  (^)  x£.  (3.20) 

In  exploring  how  to  use  this  relation,  it  is  helpful  to  note 
its  implications  in  seroth  order:  if 

ft  0  =  <>+0(t) 

n,  =  n(,0)+o(£) 

then,  since  =  0(c) 

n<0)Vx^  =  -JvxPi°0). 

This  implies  to  seroth  order,  and  conse¬ 

quently  n(,0>  =  fl(,o,(0)  as  well.  In  fact, 

n‘0,W  =  -ipl0W  (3.21) 

can  be  taken  as  a  lowest-order  characteristic  feature  of 
an  embedded  vortical  region.  It  implies  that  the  “total 
pressure*,  Px0,  tends  to  vary  only  across  (instantaneous) 
streamlines,  but  is  constant  along  them,  corresponding  to 
the  Bernoulli  result  in  lowest  order.  In  this  order,  one  could 
also  determine  by  solving 

=  -0(,0)(V>)  (3.21a) 


or 

%^  +  n,Vx^  =  -iVxPx0  (3.17) 

Ol  p 

through  first  order  in  c.  Here,  t»x  is  the  magnitude  of  £x, 
and  Px„  =  p  +  p  is  the  cross-flow  “total*  pressure  in 

the  incompressible  limit,  through  first  order  in  c. 

When  the  flow  is  irrotational,  as  in  the  "outer*  flows 
discussed  in  $  1,  Eq.  (3.17)  can  be  treated  in  the  usual 
way.  Replacing  j^x  by  Vx^  and  noting  (1,  =  0,  we  have 

<3-“> 


or 


fto  +  9±  =  Pxo(oo) 
p  dt  p 


const 


(3.19) 


for  irrotational  flow.  Note  that  in  this  frame  of  reference, 
therefore,  the  quantity  Px0  is  discontinuous  wherever  ^  or 
is  discontinuous,  as  on  the  umbilical  vortex  sheets  in  the 
wake*  and  on  the  wing  trace.  Application  of  Eq.  (3.19)  on 
both  side*  of  the  umbilical  at  the  juncture  (Fig.  6)  yields 


Aft o  _  g,  ~  ftp  _ 


dTcore 

dt 


(3.19a) 


with  appropriate  boundary  or  matching  conditions,  pro¬ 
vided  one  knew  Q(,0)(VO.  (See  [16]  and  [17]  and  5  4). 

The  presence  of  the  “perturbation*  term  involving  ££  in 
Eq.  (3.20),  however,  can  have  a  profound  influence  on  this 
conclusion,  especially  for  recirculating  embedded  flows.  For 
such  flows  the  dominant  part  of  Vx  (^).  particularly  ifO, 
is  changing  locally  in  time,  is  itself  in  the  Vxi>-direction. 
This,  in  turn,  implies  first-order  variation  of  Px0  along 
“streamlines*.  But  if  the  flow  is  in  fact  recirculating,  Px0 
must  be  cyclic  (see  §  4),  requiring  either  that  ^  =  0  in 
the  embedded  region,  or  that  an  appropriate  cyclic  average 
of  ^2-  around  the  recirculation  path  must  vanish.  Phys¬ 
ically,  for  recirculating  embedded  regions  large  enough  to 
be  essentially  inviscid,  this  observation  correspond*  to  the 
statement  that  the  local  angular  momentum  cannot  change 
other  than  by  transport  or  storage  of  that  quantity  in  a 
given  control  volume.  (Recall  that  line*  of  constant  \p,  as 
treated  hers,  are  "instantaneous  streamlines*  only:  through 
first  order  there  is  generally  a  flux  across  them.) 

To  illustrate  in  a  simple  example,  we  consider  a  circular 
embedded  region  with  solid-body  flow,  vr  =  0,  t>#  = 

Q,  =  w0(f).  where  r  is  measured  from  the  epicenter  of  the 
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convecting  core.  In  that  case,  ^  =  —s^~  +  lfo(0> 
Vx&  a  -  (*f)  t ,  where  <«,  a=  Since  <  f  x  £,  = 
-£*,  we  recover  from  Eq.  (3.30) 


U  P  _  Wo'.  .  w0r. 

-VxPio  =  — — +  —  it  ■ 


In  thie  caee,  for  an  embedded,  recirculating  circular  flow 
we  muet  conclude  that  u0  =  0,  since  =  0.  For  »uch 

a  simple  embedded  flow,  therefore,  (see  §  5)  we  would  find 
Q,  =  constant,  independent  of  time  and  space  throughout 
the  confined  region,  even  though  that  confined  region  itself 
may  grow  in  time. 


In  the  general  case,  application  of  the  cyclic  constraints 
ideas  discussed  in  $  4  provides  information  on  the  cross- 
streamline  variation  of  Pj_0 .  For  the  circular-flow  case,  as 
in  Eq.  (3.22),  the  radial  variation  of  P±„  is  determined, 
provided  we  can  determine  uo.  In  either  case,  on  matching 
the  inner  with  the  outer  flow,  one  determines  the  static 
pressure  within  the  core  as  well. 


fli0)(^)  in  (5  3.2).)  The  and  0(c)  terms  can  be 

accounted  for  by  writing 

n  =  +  (4.3) 

where  t  is  a  measure  of  the  distance  along  a  streamline,  and 
/  is  of  higher  order  than  F.  With  this  substitution  for  0,, 
Eq.  (4.2)  becomes 

^+JU.V/  =  ,V  IF;  (4.4) 

correct  to  o(^)  and  0(c).  Eq.  (4.4)  can  be  rewritten  as 

F-F,9-±.  (4.5) 

When  the  streamlines  (lines  of  constant  il>)  ‘connect*  to 
infinity,  or  to  any  region  of  known  flow  conditions,  Pj.0  (*) 
and  hence  Qj(t(>)  (=  P(^))  are  determined.  Solution  of  Eq. 
(3.13)  with  appropriate  boundary  conditions  then  proceeds 
as  usual  and  determines  the  flow  to  lowest  order. 


4  Embedded  Cross  -  Flow  Regions: 
Cyclic  Constraints,  and  The  Sheath 
Model 

4.1  Cyclic  Constraint:  Low-Shear-Stress 
Recirculating  Regions 

Determining  the  distribution  of  the  vorticity  and  the  cor¬ 
responding  total  pressure  loss  within  embedded  flow  re¬ 
quires  considerations  special  to  the  recirculatory  nature  of 
the  flows  within  them  [16],  [17].  Introduction  of  ‘cyclic 
constraints*  at  finite  Reynolds  number  to  account  for  the 
recirculating  feature  of  the  streamlines  in  the  presence  of 
viscosity  is  a  key  to  determining  D*(^)  within  the  cores. 
The  following  illustrates  the  construction  of  typical  cyclic 
constraint  conditions. 

The  full  Helmholts  equation  which  includes  viscous  ef¬ 
fects  (sometimes  called  the  dissipative  Helmholts  equation) 
is  given  by 

*Q.  =  Q,  v*,  +  ,van  (4.i) 

Using  the  slenderness  approximation,  the  z-component  of 
Eq.  (4.1)  becomes 

^  +  sivinJ  =  ^in1  +  o((1)  (4.2) 

The  first  two  terms  on  the  left-hand  side  of  (4.2)  are  of 
order  t  and  unity  respectively;  and  the  dissipation  term 
on  the  right-hand  side  is  of  0  ( ,  when  Re  being  the 
Reynolds  number.  One  solution  to  Eq.  (4.2)  is  Q,  =  0, 
representing  potential  Sow.  When  the  flow  is  rotational  and 
in  the  absence  of  strong  shear  stress,  O,  =  F(i>)  satisfies 
Eq.  (4.2)  to  seroth  order.  (Note  F(i>)  is  equivalent  to 


On  the  other  hand,  for  an  embedded  region  of  distributed 
vorticity,  the  streamlines  are  not  connected  to  an  outer, 
known  region,  and  the  above  standard  treatment  is  not  ap¬ 
plicable.  In  such  regions,  since  the  instantaneous  stream¬ 
lines  close  on  themselves  to  lowest  order,  a  different  re¬ 
quirement  must  be  met  which  will  turn  out  to  determine 
the  flow.  This  new  requirement  states  that,  for  purposes  of 
physical  continuity,  the  integral  over  a  closed  streamline  of 
the  variation  of  any  property  along  that  particular  stream¬ 
line  must  vanish.  In  other  words,  the  embedded  flow  has  to 
be  cyclic.  When  this  cyclic  constraint  is  applied  to  f{i>,  l), 
it  implies 

where  the  integral  is  to  be  taken  along  closed  streamlines. 
Eqs.  (4.5)  and  (4.6)  then  yield 


Defining  d»t  =  df/|Vj.f|,  r(tf>)  =  f  v±d»t,  r(V>)  = 
/  dsi/vx  and  also  noting  =  F'V^tp  +  E"|VxV’|a  = 
— FF'  +  v\F",  Eq.  (4.7)  becomes 

«/  [rE"  -  riff]  -  F‘  ^ j  =  0  (4-8) 

in  the  low-shear-stress  region  of  the  embedded  flow  in  the 
cross-flow  plane. 

Interpretation  of  Eq.  (4.8)  through  first  order  in  «,  i.e., 
including  the  term,  provides  additional  information. 
This  term  is  equivalent  to  the  Q  term  discussed  in  con¬ 
nection  with  the  momentum  balance  in  Eq.  (3.20)  for 
embedded  flows,  which  we  recognised  there  as  being  as¬ 
sociated  with  the  flux  of  angular  momentum  across  the 
instantaneous  recirculating  streamlines,  i>  =  const.  It  can 
be  argued,  in  fact,  that  in  large  Reynolds  number  (low- 
shear-stress)  regions  the  integrated  %  term  of  Eq.  (4.8), 


being  0(c),  if  inherently  much  larger  than  the  term  propor¬ 
tion  to  V*,  which  if  0 (]£;).  In  that  cue  we  muft  impose 
the  separate  conftnint 


Unlese  F'(i>)  if  strictly  iero  thif  implies,  a s  in  $  3,  that 
any  part  of  which  if  not  timply  periodic  in  the  cyclic 
variable  ti  muat  vanish.  Again,  for  simple  core  shapes, 
one  concludes  that  the  mean  vorticity  density  tends  to  be 
time  independent  in  embedded  recirculatory  regions.  Fi¬ 
nally,  note  that  when  Eq.  (4.8a)  applies  separately,  the 
remaining  cyclic  constraint  is  in  the  form  originally  treated 
by  Batchelor: 

i-  [r  F"  -  r  FF")  =  0  ,  (4.8b) 

corresponding  to  a  strictly  two-dimensional  embedded  vor¬ 
tical  flow  with  recirculating  streamlines. 

For  the  simple  core  geometries  that  tend  to  be  observed 
the  practical  solution  of  Eq.  (4.8b)  is  essentially  F  =  con¬ 
stant,  [16],  [17],  Only  very  weak  variation  about  this  condi¬ 
tion  can  be  found  for  reasonable  distributions  of  T(4>)  and 
r(0)  in  such  geometries.  However,  the  conclusion  F  =  con¬ 
stant,  or  essentially  so,  does  not  mean  that  the  core  flow 
is  circular  or  “solid-body*.  In  [17],  for  example,  a  range  of 
elliptical  embedded  vortical  regions,  together  with  detailed 
variations  of  F  about  its  dominant  constant  value,  were 
studied.  But  the  primary  feature  of  these  embedded  re¬ 
gions  remains  that  F  is  essentially  constant  in  space  in  the 
low-shear-stress  regime.  Our  conclusions  from  Eq.  (4.8a) 
then  provide  the  additional  information  that  F  is  constant 
in  time  as  well. 

4.2  The  Sheath  Model  and  Cyclic  Con¬ 
straints  in  General 

In  the  absence  of  any  viscosity-dominated  regions  what¬ 
soever,  Kelvin’s  theorem  would  not  allow  a  vortical  but 
essentially  inviscid  core  to  be  embedded  in  an  irrotational 
flow.  This  then  implies,  under  the  high  Re  assumption,  the 
presence  of  a  viscous  sheath  surrounding  the  core  in  order 
to  impart  vorticity  to  the  flow  entrained  into  the  core. 

Using  this  sheath  model  we  can  complete  the  determi¬ 
nation  of  the  total  pressure  distribution  in  the  cross-flow 
plane.  On  account  of  Eqs.  (3.19)  and  (3.20),  a  difference 
in  the  evolution  of  total  pressure,  P±0,  inside  and  outside 
the  core  is  inevitable.  The  sheath  plays  the  role  of  redis¬ 
tributing  Pjl0  at  the  core  edges,  thus  smoothing  out  any 
otherwise  forbidden  “jumps*  in  Pj.0-  Moreover,  the  length 
around  the  perimeter  of  the  cores  needed  for  the  viscous  ac¬ 
tion  to  provide  a  specified  change  of  P±0  grows  with  some 
power  of  the  Reynolds  number  [17].  This  fact,  combined 
with  the  required  sheet  strengths  of  the  umbilicals,  will  be 
used  in  $  5  to  determine  the  core  sites  as  they  develop. 

An  example  involving  both  the  application  of  Batche¬ 
lor’s  cyclic  constraint  and  the  redistribution  of  total  pres¬ 
sure  by  the  viscous  sheaths  is  illustrated  in  (17|.  In  that 


work  a  model  for  the  leading  edge  separation  and  reattach¬ 
ment  over  a  two-dimensional  plate  is  developed.  The  study 
utilises  Eq.  (4.8b)  to  deduce  the  form  of  vorticity  in  the 
embedded  core  regions.  The  model  for  the  viscous  sheath 
consists  of  a  shear  layer  made  up  of  two  sub-layers.  (See 
Fig.  7)  The  inner  sublayer,  being  in  direct  contact  with  the 
embedded  region,  consists  of  recirculating  flow  that  incor¬ 
porates  the  boundary  layer  on  the  upper-side  of  the  plate. 
In  contrast  the  outer  sub-layer,  of  which  a  portion  is  the 
boundary  layer  on  the  underside  of  the  plate,  is  bounded 
by  the  irrotational  outer  flow. 

In  this  sheath  model  the  outer  sub-layer  is  considered 
to  drive  the  inner  sub-layer  via  shear  stress.  One  impor¬ 
tant  criterion  emerging  from  this  two-layer  model  is  that,  to 
maintain  its  recirculating  nature,  the  inner  sub-layer  flow 
must  be  given  enough  energy  to  overcome  viscous  dissipa¬ 
tion  on  the  plate.  In  other  words,  the  gain  in  total  pressure 
of  a  fluid  particle  travelling  from  A  through  B  to  C  must 
turn  out  to  exactly  balance  the  loes  in  the  same  quantity 
as  the  particle  completes  the  cycle  by  moving  from  C  back 
to  A.  This  corresponds  to  an  application  of  the  general¬ 
ised  cyclic  constraint  idea.  It  is  used  together  with  the  fact 
that  the  static  pressure  along  the  outer  edge  of  the  sheath 
is  governed  by  the  outer  flow  to  deduce  the  sheath  profile 
and  its  peripheral  static  pressure  distribution  via  an  iter¬ 
ative  procedure.  Figs.  8  and  9  show  results  obtained  for 
two  different  Reynolds  numbers:  streamlines  for  the  exter¬ 
nal  flow  are  presented  in  (a)  and  those  for  the  embedded 
region  in  (b).  Such  results  provide  a  framework  from  which 
techniques  developed  so  far  can  be  extended  to  include  cor¬ 
rections  for  unsteady  effects.  But  in  addition,  as  illustrated 
in  5  5,  these  ideas  provide  one  of  the  Anal  links  necessary 
to  determine  the  vortex  core  evolution. 

5  Closure  of  the  Aerodynamic  Prob¬ 
lem  with  an  Embedded  Circular 
Core 

This  section  illustrates  how  one  can  systematically  em¬ 
ploy  the  techniques  developed  in  the  previous  sections  to 
determine  the  flow  in  an  embedded  region  based  on  a  sim¬ 
ple,  but  physically  revealing  model.  In  particular  the  role 
of  interaction  between  the  umbilical  and  viscous  sheath  is 
clearly  demonstrated.  The  model  consists  of  a  strictly  cir¬ 
cular  core,  with  a  convective  epicenter  and  with  a  spatially 
constant,  time-invariant  vorticity  density  w0-  This  vorticity 
distribution  is  consistent  with  the  results  of  5  4.  A  “core* 
can  be  said  to  be  present  when  the  structure  of  a  portion 
of  the  wake  becomes  so  rolled  up  as  to  form  a  spiral  with 
spacing  between  successive  turns  which  is  comparable  to 
that  of  an  appropriate  diffusion  thickness  of  the  umbilical 
sheet  feeding  it. 

Using  a  simple  estimate  for  the  influence  of  viscosity  on 
the  umbilical  sheet  we  find  that  each  circulation  element 
has  a  thickness  which  develops  as 

6  =  4  \/ut 


(5.1) 


where  v  i*  the  kinematic  viscosity  of  the  fluid  and  r  i*  the 
'drift  time*  aeeociated  with  that  particular  v-ake  element. 
It  follow*  from  (S.l)  that  6  lucre  a*  ee  monotonically  along 
the  umbilical  from  *«ro  thickneee  at  the  wing  edge,  to  a 
maximum  at  the  point  where  the  umbilical  join*  the  vi*cou* 
theath  surrounding  the  core.  The  drift  time  at  thi*  point 
i*  ’junct  =  *  "  ^junct  (5  2).  Alternatively  (5.1)  can  be 


where  R,c  i*  the  Reynold*  number  baeed  on  the  wing  chord 
length,  R,c  =  t/*a-c|/°* a* ,  and  To  i*  the  characterise  time 
of  the  un*teady  crossflow  problem  T0  =  a  • 

Having  introduced  the  concept  of  a  diffusion  thickne**, 
we  proceed  to  eetabliah  a  “cut-off*  criterion  for  determin¬ 
ing  the  core  formation.  We  firtt  approximate  the  rolled-up 
portion  of  the  ideal  wake  by  a  tpiral  defined  by 

r(9)  =  ~  (5.2) 

where  r(i)  i*  the  diatance  from  the  center  of  the  apiral  at 
angle  8,  and  r0  i*  a  conatant.  Denoting  the  number  of 
encirclement*  about  the  epiral  centre  by  N,  the  relation 

r»  =  N  to  (5.3) 

ia  obtained  for  N  >>  1.  We  can  now  seek  a  relation  be¬ 
tween  r0  and  £junct  with  the  aim  of  establishing  the  cut-off 
criterion  for  the  core;  that  is,  determining  where  and  when 
to  replace  the  ideal  wake  by  an  equivalent  vortical  core.  We 
choose 

r°  ~  ^junct  (®-*) 

since  in  that  case  one  can  no  longer  ignore  the  effect*  re¬ 
sulting  from  the  diffusion  of  vorticity  among  the  vortex  el¬ 
ements  in  the  contorted  wake  region,  i.e.,  the  core  becomes 
smeared  out  and  diffused. 

Next,  recall  from  §  3  that  once  a  circular  core  is  estab¬ 
lished  at  tome  initial  stage,  it  will  have  a  spatially-constant, 
time-invariant  vorticity,  wo.  The  drift  time  can  now  be  em¬ 
ployed  to  formulate 

form)  d8 

&T-t  fyorm  Hunct  ~  t  u\ 


/•’“("-"form)  J8__ 

"  Jo  m 
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junct  (*  ~  ‘form  + 
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Additional  information  relating  »junct,  "o  and  r»  can  in 
principle  be  obtained  from  a  sheath  model  like  that  of  Ref. 
[17] .  However  in  this  context  it  suffices  to  adopt  an  even 
simpler  model  in  which  the  theath  is  taken  to  be  of  constant 
thickness  equal  to  that  of  the  umbilical  at  ’)uncture.  In  this 
case  the  total  pressure  lose  around  the  perimeter  of  the  core 
is  given  by 

-Apo=^-Af  (5.7) 

fi  junct 

where  U,  is  the  velocity  at  the  edge  of  the  core  U.  = 
and  At  is  the  perimeter  of  the  core,  At  =  2 xr».  Recalling 
from  Section  (3.2)  that  the  total  pressure  lose  around  the 
core  is  equal  to  the  rate  of  change  of  core  circulation,  and 
that  for  a  solid  body  core, 


Tcore  =  x  w0  rj  , 


we  finally  obtain 


where  tform  is  the  time  when  the  core  first  appears;  /Vforal 
corresponds  to  ff  at  t  =  fform.  Recall  that  nct  i*  the 
drift  time  associated  with  the  vortex  element  at  the  umbil¬ 
ical/core  juncture  at  time  (.  Eqs.  (5.3),  (5.4),  (5.5)  can  be 
combined  with  To  and  R,c  to  give 


-(?) 


In  the  long  run,  t  »  (-tfonn  +  jftyorm)-  (5-6) 
becomes 

the  unknown  parameters  being  Tjunct,  rs  and  w0. 

Inspection  of  Eq.  (5.10)  reveals  two  limits  of  the  behavior 
of  tjUnct  which  give  important  insight  into  the  rate  of  core 
growth.  The  first  is  the  one  for  which 


This  is  the  limit  of  a  small  or  very  slowly  growing  core.  In 
this  case,  using  (5.9)  the  core  grows  asymptotically  as 

r* - —t*  .  (5.12) 

V 

The  second  limit  illustrates  the  case  in  which  the  core  is 
entraining  wake  vorticity  as  fast  as  the  umbilical  can  supply 
it.  For  sufficiently  large  times, 


and  the  core  grows  according  to 

r‘ - 7==*^  ■  (51<) 

Eq.  (2.16)  allows  one  to  express  at  any  instant  the  total 
circulation  of  the  vortex  elements  lying  between  any  loca¬ 
tion  along  the  wake  and  the  free  end  of  the  wake  a*  a  func¬ 
tion  of  the  drift  time  associated  with  the  element  at  that 
particular  location.  Hence  there  exist*  a  one-to-one  corre¬ 
spondence  between  Tcore  and  Tjanct.  Eliminating  wo  from 
Eq*.  (5.8)  and  (5.10)  then  gives  qua  function  of  fjunct- 
Such  a  relation  between  r»  and  rjunct,  together  with  Eq. 
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(S.9),  finally  determine*  ^unc*  and  r».  A*  a  re*ult  Tcore 
and  uo  are  determined. 

The  analyti*  can  al*o  be  extended  to  determine  Njorm 
aince  wo  i*  now  known  and  conataat  in  time.  Thi*  it 
achieved  by  projecting  backwards  in  time  to  f  =  *form' 
Eqs.  (S.6),  (5.8),  (5.9)  are  then  solved  in  much  the  same 
way  as  before,  the  anknowns  now  being  N{orm,  r»,  rjunc{ 
instead  of  w0,  r»,  ^uact.  Once  lVform  is  obtained,  t{orm  is 
determined. 

Once  the  shape,  location  and  circulation  distribution  of 
each  wake  has  been  determined,  the  information  is  at  hand 
to  calculate  the  loading  and  force  on  the  wing.  Eq.  (2.15) 
is  used  to  find  the  induced  velocity  at  the  wing  trace  associ¬ 
ated  with  the  wake.  Eq.  (2.9)  then  serves  to  determine  the 
discontinuous  part  of  the  tangential  velocity  required  to 
satisfy  the  boundary  condition  Eq.  (2.7)  of  no  flow  through 
the  trace.  This  is  added  to  the  ‘‘wake  free*  velocity  jump 
AVo  given  by  Eq.  (2.8).  Then  the  unsteady  Bernoulli  Equa¬ 
tion  determine*  the  instantaneous  loading  at  each  point  on 
the  trace,  (as  in  |6j)  which  can  be  integrated  to  yield  the 
normal  force  on  the  wing  per  unit  chord.  Hence  the  local 
force  coefficient  (i.e.,  the  force  per  unit  chord  normalised 
by  the  dynamic  pressure  and  the  local  span)  is  given  by 


6.  McCune,  J.E.,  “Interactive  Aerodynamics  of  Wings 
in  Severe  Maneuver."  To  be  included  in  Proceeding* 
of  Workshop  D  on  Unsteady  Flow,  United  States  Air 
Force  Academy,  July  28-30,  1987.  To  be  published. 

7.  Hoeijmakers,  H.W.M.,  Numerical  Simulation  of  Vorti¬ 
cal  Flow,  MP  86032  U,  NLR,  1986. 

8.  Powell,  K.G.,  "Vortical  Solutions  of  the  Conical  Euler 
Equations."  M.I.T.  Ph.D.  Thesis,  June  1987. 

9.  Earnshaw,  P.B.,  "An  Experimental  Investigation  of  the 
Structure  of  a  Leading  Edge  Vortex."  Aero.  Res. 
Coun.,  Lond.,  Rep.  No.  22,  876. 

10.  Thomas,  J.L.,  Taylor,  S.L.  Anderson,  W.K..  "Navier- 
Stokes  Computations  of  Vortical  Flows  Over  Low  As¬ 
pect  Ratio  Wings."  AiAA  87-0207,  January  1987. 

11.  Newsome,  R.W.,  Kandil,  O.A.,  "Vortical  Flow  Aero¬ 
dynamics  —  Physical  Aspects  and  Numerical  Simula¬ 
tion."  AIAA  87-0205,  January  1987. 

12.  Munk,  M.,  "The  Aerodynamic  Forces  on  Airship 
Hulls."  NACA  Report  184. 


v/4? 

Cn  =  — —  sin  a  +  wake  history  effects 


(5.15) 


where  the  first  term  is  Jones’  classical  wake  free  contribu¬ 
tion  and  where  the  wake  history  effects  are  calculated  by 
the  integrals  just  cited. 

This  result  emphasises  that  the  forces  depend  not  only 
on  the  local  angle  of  attack  and  spanwise  growth  of  the 
trace  but  also  on  the  cumulative  effects  of  all  events  in  the 
evolution  of  the  wakes  up  to  that  point.  This  would  imply 

change  in  the  magnitude  and  distribution  of  the  normal 
force  across  the  chord  which  will  have  a  substantial  effect 
on  the  wing’s  overall  lift,  drag,  and  moment. 
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